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Abstract 

This contribution reviews the present status of the Skyrme-Hartree- 
Fock (SHF) approach as one of the leading self-consistent mean-field 
models in the physics of atomic nuclei. It starts with a brief sum¬ 
mary of the formalism and strategy for proper calibration of the 
SHF functional. The main emphasis lies on an exploration of the re¬ 
liability of predictions, particularly in the regime of extrapolations. 
Various strategies are discussed to explore the statistical and system¬ 
atic errors of SHF. The strategies are illustrated on examples from 
actual applications. Variations of model and fit data are used to get 
an idea about systematic errors. The statistical error is evaluated in 
straightforward manner by statistical analysis based on fits. This 
also allows also to evaluate the correlations (covariances) between 
observables which provides useful insights into the structure of the 
model and of the fitting strategy. 

1. Introduction 

This contribution is devoted to a brief review of the 
Skyrme-Hartree-Fock (SHF) approach to nuclear struc¬ 
ture and dynamics with emphasis on a critical analysis of 
its predictive power. SHF belongs to the class of nuclear 
self-consistent mean-field models which became competi¬ 
tive in the early 1970ies. The first nuclear models were 
purely macroscopic of the type of a liquid drop model. 
The observation of shell structure led to the invention of 
an empirical mean-field model, the nuclear shell model, 
with which one could reveal the basic mechanisms of nu¬ 
clear structure from magic nuclei over deformation to low- 
energy collective motion for overviews see the much cele¬ 
brated nobel lectures mw and books [4] . These findings 
have confirmed that nuclear structure and dynamics can 
be well described in terms of a mean-field. The mixed 
macroscopic-microscopic model employing empirical shell 
models has been driven to a precise instrument of describ¬ 
ing and predicting nuclear masses, see e.g. [5]. Paral¬ 
lel to these developments, one has looked for the natural 
next step, namely for self-consistent mean-field models of 
the type of the Hartree-Fock approximation. Self consis¬ 
tency eliminates the need for intuitive explicit shaping of 
the mean-field potential and is thus applicable to a much 
broader range of processes. 

The SHF approach resembles Hartree-Fock but uses an 
effective interaction, called the Skyrme force. In fact, it is 
more precise to say that it uses an effective energy func¬ 
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tional derived from the Skyrme force. The Skyrme force 
was proposed already in [5] , had a break-through for prac¬ 
tical applications with the first high-quality parametriza- 
tions coming up [7], and has evolved meanwhile to one of 
the most widely used standard schemes for self-consistent 
nuclear modeling. There are many more models of that 
sort which rely on effective energy functionals. The two 
strongest competitors of SHF are the Gogny force UM 
and the relativistic mean-field (RMF) model [101 El] for 
which meanwhile a couple of variants exists (TO). For a 
review of all three nuclear effective energy functionals see 
[TO] . It is interesting to note that these nuclear effective 
energy functionals became fashionable almost at the same 
time as density functional theory for electronic systems 
[HEilEl]- The difference is that electronic energy func¬ 
tionals have been derived from first principle and ab-initio 
calculations while nuclear functionals are mostly deter¬ 
mined by calibration with respect to empirical data. That 
is because the nuclear many-body problem is by orders 
of magnitude more involved than that for electrons. In 
spite of considerable progress in nuclear ab-initio calcula¬ 
tions (for reviews see e.g. [TTl EHl EH EQ] E]), the case 
is still not sufficiently well settled to serve as basis for 
high-quality functionals. The standard pathway is since 
the early days to motivate the form of the functional by 
the basic structure of a low-momentum expansion and to 
calibrate the parameters of the functional by fits to em¬ 
pirical data, usually nuclear ground state properties. This 
raises the problem how to get an idea about the predic¬ 
tive power of the effective energy functional and to es¬ 
timate the uncertainties in extrapolations to observables 
outside the pool of fit data. As we know, it holds for every 
physical theory that it can never be ultimately proven, but 
only be confirmed stepwise by experience until one hits the 
range of validity. This holds even more true for nuclear 
mean-field models in their pragmatic mix of theoretical 
and empirical input. It is the main topic of this contri¬ 
bution to present and discuss the various approaches to 
estimating the uncertainties within nuclear self-consistent 
mean-field models. Thereby we confine the considerations 
to the SHF functional as a prototype example. The results 
are transferable to the competitors, the Gogny force and 
the RMF. 

The paper is outlined as follows: Section [^ summarizes 
the SHF functional as it is used throughout the paper 
and sketches its formal motivation. Section [^ briefly col¬ 
lects the various observables employed for calibration of 
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the SHF functional (bulk properties of the ground states 
of finite nuclei), for characterization (response properties 
of homogeneous nuclear matter), and in a broad range of 
further applications. Section introduces the method of 
least-squares (x^) fits and the subsequent statistical analy¬ 
sis delivering extrapolation uncertainties and correlations 
(also coined covariances) between observables. Section]^ 
addresses the problem of error estimates from many differ¬ 
ent perspectives. Section exemplifies most of the meth¬ 
ods compiled in sectionj^in terms of practical applications 
to a variety of different observables. 


2. The SHF energy-density functional 

2.1. Basic constituents: Densities and currents 

Mean-field theories describe a system through a set of 
single particle (s.p.) wavefunctions ifiair) and associated 
BCS amplitudes Va for occupation and the complementing 
Ua = \/l — r’a for non-occupation. These are summarized 
in the one-body density matrix 

Qqi^y) ='^Wavlip^{Y)ip\{r') (1) 

a^q 

where q labels the nucleon species with q = p for protons 
and q = n for neutrons. The Wa is a further factor which 
serves to terminate the summations, so to say a cutoff for 
pairing space. It will be discussed in connection with the 
pairing functional in section [2. 2. 3[ 

The Skyrme-Hartree-Fock (SHF) energy-density func¬ 
tional requires the knowledge of only a few local densities 
and and currents. These are, sorted according to time 
parity: 


time even: 

Pq = tr<,{p(r,r')}|^^^, 

Tq = tr<,{V^Vr/p(r,r')}| , 

Jq = -itr^{VrXdp(r,r')}[^^^, 


= density 
= kinetic density 
= spin-orbit density 


time odd: 

(Tq = tr<,{<Tp(r,r')}|^^^, 

jq = 3m{tr„{V^£<(r,r')}|^^^,} 

Tq = -tr,^{dVrV^/6i(r,r')}|^^^, 


= spin density 
= current 

= kinetic spin-dens. 


time mixed: 

^q = = pairing density 

( 2 ) 

It is advantageous to handle the densities in terms of 
isospin T S {0,1} instead of protons p and neutrons n. 
Thus we will often consider the recoupled forms which 
read for the local density 


PO — P — Pp “f Pn ) Pi — Pp Pn ) (3) 


and similarly for the other densities and currents. The 
isoscalar density po = P is equivalent to the total density 
and the difference pi corresponds to the isovector density. 
All densities and currents in the collection ([^ are real and 
have definite time parity, except for the pairing density ^ 
which is complex and has mixed time parity. 


Note that the above collection of densities and cur¬ 
rents is richer than in electronic density-functional the¬ 
ory (DFT) which employs usually only the local (spin) 
density Pai^) [22]. This indicates that nuclear DFT is 
different [23) . This will also become apparent later on in 
connection with the calibration of the energy functional, 
see section which proceeds much different than for elec¬ 
tronic DFT. 


2.2. The composition of the total energy 

Starting point of all self-consistent mean-field theories 
based on DFT is an expression for the total energy. This 
reads for SHF 


^total — I d T ^'kin j d T ^Sk -^Coul ^pair 


The one-body kinetic energy is given by 

^kin — 'Yv ^ 'Yn 


2m„ 


2m„ 


(4) 


(5) 


Keeping it at the exact level implies that full quantum 
mechanical shell structure is maintained which means that 
we deal with DFT at the level of the Kohn-Sham approach 

[mnn. 


2.2.1. The Skyrme energy functional 

Key piece is the SHF energy density: 

^Sk ~ ^Sk,even T ^Sk,odd 5 

^Sk,even — 


(6a) 
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4Cf(Ti-ri 


(6c) 

The formal reasoning for this functional will be presented 
in section 2.3 But already here, the building principle is 


obvious: The energy density contains all conceivable bi¬ 
linear couplings of densities and currents up to second or¬ 
der in derivatives. Each term is time even, although it may 
be composed of a product of two time-odd currents. The 
coupling constants are denoted as with obvious ab¬ 

breviations for “type”. Each coupling constant depends, 
in principle, on the isoscalar density po- In practice, one 
takes a minimalistic approach and augments only the pp 
term by a minimalistic density dependence oc (7^’“. This 
suffices to deliver a good description of ground state prop¬ 
erties and excitations. On the other hand, more density 
dependence is hard to determine empirically from nuclear 
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data because finite nuclei cover a small range of densities 
due to nuclear saturation. 

There is a hierarchy of importance in the various terms 
of the functional Q. The leading terms are those which 
depend on px only. These already provide a good descrip¬ 
tion of bulk matter and an acceptable, although rough, 
zeroth order description of nuclei, see section |6.1| Terms 
with derivatives add details which lift the model to the 
level of a quantitative description. The minimal set of 
terms is indicated by yellow (gray) shading. Already these 
terms alone provide an excellent description of bulk prop¬ 
erties of hnite nuclei (energy, charge radius, charge surface 
thickness). All further terms are necessary to allow also 
a good description of response properties as giant reso¬ 
nances, polarizability, low-energy vibrations, or fission. A 
large part of the discussions in this manuscript deals with 
exploring the impact of the various terms on the modeling. 


2.2.2. Coulomb energy 

The Coulomb energy is treated as 

.^Coul — .^Coul.dir .^Coul.ex ; (^^) 



where = 1.44MeVfm. Coulomb exchange is done at 
the level of the local-density approximation [52]. It is to 
be noted that not all mean-field models include Coulomb 
exchange (for example, most RMF models do not), in the 
present examples, it is included. The direct Coulomb en¬ 
ergy should employ, in principle, the charge density pc- 
However, the mere proton density Pp is used in most ap¬ 
plications. This makes not much of a difference and is 
simpler to handle. We follow this tradition also here. 


The part fsk.odd collecting all coupling with time-odd 
currents is inactive in static calculations, e.g., for ground 
states of even-even nuclei. They come into play with exci¬ 
tations and with odd nuclei. Note that part of the coeffi¬ 
cients, those printed in gray, are taken over from the time- 
even part. There is no freedom to chose them differently 
because identity of these coefficients is crucial to guarantee 
Galilean invariance of the functional [ISIEl]. There are 
choices for the other coefficients. They are all connected 
with terms carrying spin and become active only only for 
odd nuclei and spin excitations. We will postpone a dis¬ 
cussion of these spin terms and their coefficients to section 
16.61 All excitation modes in even-even nuclei with natu¬ 
ral parity (giant resonances, low-lying vibrational states, 
rotation) access only the terms in £sk,odd which are fixed 
by Galilean invariance. 


The functional (6b I contains two terms leading to a 
spin-orbit potential. The compulsory basic term is the 
one oc pV-J. Its strength is an independent param¬ 
eter of the SHF functional (much unlike relativistic mod¬ 
els where the spin-orbit strength is automatically given 
by the Dirac structure of the wavefunctions uni dn ESI). 
The other term oc is called tensor spin-orbit term. Its 
parameter is linked to the parameter of the kinetic 
term, see table [^ In fact, we quote here only a simplified 
version. The full tensor term has true tensor structure 
(not reduced to a vector J) and there is a complementary 
tensor term. These tensor terms are ignored in the major¬ 
ity of SHF applications, i.e. one deals with c;f = 0. This 
is what we will assume throughout the following. We refer 
to [5H][57] for a very detailed survey of tensor terms and 
its effects in a great variety of observables. 


Before continuing with the other contributions to the 
total energy Q , we mention in passing that the SHF func¬ 
tional ^ does not necessarily guarantee unconditional 
stability. In particular spin and tensor terms are likely 
to induce instabilities in symmetry unrestricted calcula¬ 
tions. This important issue has been discussed over the 
years in several respects. For recent extensive studies see, 
e.g., [Ml EH EH ESI . 


2.2.3. Pairing energy 


Pairing is inevitable in case of open shell nuclei. They 
have a high density of almost degenerated states. This 
gives the residual two-body interaction a chance to mix 
these states in order to produce a unique ground state 
[30] . Pairing explores what is called the particle-particle 
(pp) channel of the effective interaction which is differ¬ 
ent from the particle-hole {ph) channel responsible for the 
mean field m- Thus it is justified and customary to use 
a separate functional for pairing, namely 


Apair — . 


E 

qe{p,n} 




pair,q 





Po 

POjpair 


( 8 ) 


where p and ^ are the particle and pairing densities as 
defined in eq. (HI- The pairing functional (§ contains a 
continuous switch, the parameter po.pair, where a pure <5- 
interaction (DI) is recovered for po.pair —> oo also called 
volume pairing. The general case is the density depen¬ 
dent d-interaction (DDDI). A typical value near matter 
equilibrium density po.pair = 0.16 fm“^ concentrates pair¬ 
ing to the surface. This is often denoted as surface pair¬ 
ing. Allowing po.pair to be a free parameter of the model 
yields a value po.pair ~ 0.2 fm“^ which puts the actual 
functional somewhere half way between volume and sur¬ 
face pairing [32]. A fully variational treatment of pairing 
yields the Hartree-Fock-Bogoliubov (HFB) equations. A 
widely used and robust simplification is the BCS approx¬ 
imation which suffices for all well bound nuclei [33]. All 
following calculations are done within the BCS scheme. 

The space of pairing-active states has to be limited be¬ 
cause the zero-range nature of the pairing functional (|^ 
produces poor convergence with size of phase space. In nu¬ 
merical calculations only a moderately large set of states 
can be included. This is expressed by the phase-space 
weight Wa in the dehnition Q of the on-body density. 
Older recipes employ a sharp cutoff Wa{sa) = ff(ecut — £a)- 
This can be done in connection with large pairing spaces 
reaching up to 50 MeV above the Fermi surface [34] . 
Smaller spaces require a smooth cut off for which one often 
uses 


Wa = [1 -b exp ((£a - {ep + ecut))/Ae)] ^ (9) 
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where typically Ccut = 5MeV and Ae = Ccut/lO [551155] . 
The value Cq is the single particle energy of the state a 
and Ci? is the chemical potential. This works very well 
for all stable and moderately exotic nuclei. For better 
extrapolation ability away from the valley of stability, the 
fixed margin Ccut is modified to use a band of fixed particle 
number oc instead of a fixed energy band m- 

2.2.4. Correlation energy 

Finally, we come to the correlation energy. It contains 
several contributions 

^corr — F^cm T ^rot T .^vib ; 

a correction Ecm for the spurious center-of-mass energy, 
another term i?rot for rotational projection (non-zero only 
in deformed nuclei), and a term E^ih accounting for soft 
surface vibrations. The leading part which is used in prac¬ 
tically all applications is i?cm- There exist several variants 
for it as summarized, e.g., in miMiisn]. We use it here 
in the form 

Ecm = 2mA ’ Pn (Ha) 

n—1 

which accounts for the effect of center-of-mass projection 
evaluated in second order Gaussian overlap approximation 
[40] . It employs, in fact, a two-body operator which makes 
its application in variational equations very cumbersome. 
It is thus only evaluated a posteriori, i.e. for the given 
solution of the SHF mean-field equations. 

The rotational correction is, similar as Ecm, an approx¬ 
imation to rotational projection and looks 


and because the (heavily fluctuating) correlations from 
low-energy vibrations cannot be embodied into a smooth 
energy functional. All three of these correlations are as¬ 
sociated with low-energy, large-amplitude collective mo¬ 
tion. There emerges now a subtle problem with extend¬ 
ing the SHF functional to these cases. It is by defini¬ 
tion a density functional, well defined only as expectation 
value over one mean-field state producing unambiguously 
one density matrix Q and local densities therefrom. But 
the generator-coordinate method of large amplitude col¬ 
lective motion requires energy overlaps between different 
mean-field states [15]. These are a priori undefined for 
the SHF functional. One can motivate a unique extension 
to compute those overlaps which works well for center-of- 
mass projection [40]. But it was figured out later that 
one runs into subtle problems with the analytical struc¬ 
ture of such an extended definition particularly in cases 
of particle-number projection ST], but also for rotational 
projection. Luckily enough, these problems are still ab¬ 
sent when computing energy overlaps of mean-field states 
which stay still close to each other. And only these near 
overlaps are employed in modeling large-amplitude collec¬ 
tive motion within the Gaussian overlap approximation to 
the generate-coordinate method [25]. All the correlations 


required here (see section 2.2.4) are evaluated in this limit 


and are thus at the safe side. The problem persists when 
full projection is necessary. One solution is to develop an 
effective energy functional which is developed consistently 
from an effective interaction. Work in this direction is in 
progress, see e.g. [55] . 


2.3. Motivation of the SHF functional 


Ernt — 


20 4 


( 12 ) 


g{x) = xda: log ( I dye 




where J is the operator of total angular momentum and 
0 the momentum of inertia as evaluated in the ATDHF 
approximation (also coined self-consistent Inglis inertia) 
[13 E]. This rotational correction plays a crucial role for 
well deformed nuclei and should be included for them. The 
switch factor g serves to limit the correction to deformed 
nuclei and to leave (nearly) spherical nuclei untouched [5T1 


Similar as the c.m. correction E, 


(full) 


it is a two- 


body operator and thus only evaluated a posteriori. 

Finally, i?vib becomes noticeable in all nuclei with soft 
surface vibrations which are typically the transitional nu¬ 
clei between spherical and well deformed ones. Its evalu¬ 
ation is very involved, for details see [ 111123125 ]. We will 
include the full correlation energies only once in the in¬ 
troductory overview in section 6.0.1] All further examples 
are confined to nuclei which are proven to have negligible 
vibrational-rotational corrections and need only Ecm- 


2.2.5. A comment on large amplitude motion 

The correlation energies discussed in the previous section 


2.2.4 are, in fact, covering effects beyond mean field. They 
are included for symmetry reasons (translation, rotation) 


The standard pathway to derive an energy-density func¬ 
tional from ab-initio calculations is for electronic systems 
the much celebrated local-density approximation (LDA) 
[55] refined by non-local corrections through the gener¬ 
alized gradient approximation (GGA) [T^- This line of 
development is also a strong motivation for the nuclear 
case. But it cannot be simply copied for nuclear systems, 
first because nuclear ab-initio calculations are not yet as 
reliable as electronic ones are, and second, because nu¬ 
clear energy functionals require more than mere density 
dependence as can be seen from the SHF functional ([^. 
Although a quantitative derivation from ab-initio theories 
is still inhibited, at least the desirable formal structure 
of a Skyrme-like effective energy functional can be moti¬ 
vated by the technique of the density-matrix expansion 
[49l [50] . We sketch here the basic steps. For simplicity, 
we concentrate on the spatial part of the expansion and 
ignore spin-isospin structure. We also skip explicit vector 
notation. 

Most ab-initio models deliver at the end an effective two- 
body interaction for an underlying mean-field calculation 
in terms of an involved integral operator, the T-matrix 
T. An example is the Brueckner-Hartree-Fock method 
(BHF) whose T-matrix serves finally as effective force for 
the Hartree-Fock part, for reviews see e.g m- An effec¬ 
tive force as integral operator can also be extracted from 
other ab-initio models m as, e.g., the unitary correlator 
method. In any case, the most general total interaction 
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energy reads 

■Spot = Jdy'eix, x') T{x, x'; y, y') g{y, y') (13) 

where g{x, x') is the one-body density matrix Q. The key 
point is that g(x^ x') varies slowly with a;, x' within a typ¬ 
ical range of where kp is the Fermi momentum. The 
T-matrix, on the other hand, is well concentrated in space, 
non-zero only for small differences in all pairs of coordi¬ 
nates with typical ranges mostly below kp^. This suggests 
a Taylor expansion in all four coordinates x, x', y, y' about 
the common center R = (x-\-x' + y + y')/A. This reads up 
to second order 


g{x,x') « p{R) + {x — R)'^p\^ +\{x — Kf/S.p 


1 / /, 9 / ^ A 1 I 

2V- - / V ^^P)\r 


+i{x-x')j\j^ + -{x-x'f{T--Ap)\i^ (14) 


where we abbreviate x = {x + x')l2. A similar expan¬ 
sion holds for g{y,y'). We insert this into the interaction 
energy (13) keep all terms up to second order in deriva¬ 


tives and recall that the T-matrix conserves parity as well 
as time-parity. This eliminates all terms of first order. 
What remains is just the SHF functional of course, 
without the spin and spin-orbit terms which had been ig¬ 
nored in this quick “derivation”. The basic ingredient, the 
T-matrix, is assumed to stem from homogeneous matter 
and thus depends on the density for which it was evalu¬ 
ated. Thus all expansion coefficients (7^^® in the ansatz 
([^ carry, in principle, a density dependence. This is as 
far as we can get with formal reasoning. As argued above, 
we are still lacking sufficient information to fix all possible 
density dependencies and associate one only to the lead¬ 
ing zeroth order terms oc Pq and oc pp This looks at first 
glance like a somewhat helpless escape. But it turns out to 
be a pragmatic and fruitful guess delivering a high-quality 
functional for many purposes, as we will see. The reason is 
that the density of nuclei gathers predominantly around 
the equilibrium density po « Peq due to strong nuclear 
saturation. 

Altogether, the density-matrix expansion demonstrates 
how a zero-range effective interaction emerges naturally 
from the initially given involved operator structure. We 
are dealing with a typical low-energy or low-g expansion 
m- It requires that the spatial structure should be suffi¬ 
ciently smooth which means in the present case that the 
essential physics of the T-matrix is concentrated at length 
scales smaller than the typical wavelength kj,^. This is 
fine at normal nuclear density. But one should warned 
with extensions to high densities. Sooner or later, effec¬ 
tive functionals of zero range will become inappropriate 
although the validity of effective functionals exceeds often 
the range of such safe estimates. 


3. Observables 

3.1. Homogeneous nuclear matter 

Although not experimentally accessible, homogeneous nu¬ 
clear matter is an extremely useful system to character¬ 
ize basic nuclear properties. Nuclear matter properties 
(NMP) have been used over decades as key parameters 
in macroscopic models [3 [521 [53]. Extensive studies in 


this domain have developed reliable, although model de¬ 
pendent, values for them. It is technically simple for ev¬ 
ery mean-field model to compute NMP for homogeneous 
matter. Moreover, NMP can be given a physical interpre¬ 
tation. Thus NMP are helpful quantities to characterize 
a model and to compare different models. We introduce 
here those NMP which will be used in the following. 

Homogeneous nuclear matter is to be taken without 
Coulomb force (it can be assumed to be neutralized by 
a homogeneous electron background), pairing, and c.m. 
correction. It remains the energy per particle as 

-r[Po,Pi,To,Ti) = - . (15) 

A Po 

where we consider for a while p and r as independent 
variables. We keep the independence just to allow a 
simple computation of effective masses. Of course, a 
given system is characterized just by the densities pp and 
the kinetic density depends on these given densities as 
tt = tt{po,Pi)- Thus we have to distinguish between 
partial derivatives d/dr which take Tp as independent and 
total derivatives d/dp which know only the pp dependence. 
The relation is 

— = —061 

dpT dpp “ dpp drpi 


The standard NMP are defined at the equilibrium point 
{po = Peq, Pi = 0) of symmetric nuclear matter. They are 
summarized in tabled) Most of the NMP are obvious. 

A few remarks on two more subtle features: The slope 
of symmetry energy L characterizes the density depen¬ 
dence of the symmetry energy which allows to estimate 
the symmetry energy at half density, i.e. at surface of 
finite nuclei. The enhancement factor for the Thomas- 
Reiche-Kuhn (TRK) sum rule [33] is a widely used way to 
characterize the isovector effective mass which is obvious 
from the given expression involving derivative with respect 


to Tl. 


The NMP in table [T] can be grouped into four classes: 
first, the (isoscalar) ground state properties Peq and 

E/A , second, isoscalar response properties K and m/m, 


and third, isovector response properties J, L, ktrk- The 
response properties determine zero sound in matter [55] 
and subsequently they are closely related to giant reso¬ 
nance modes in finite nuclei as we will see later. There is 
a fourth category, the surface energies which go already be¬ 
yond homogeneous matter and are explored in the surface 
of semi-infinite matter. Their computation in the context 
of quantum mechanical mean-field theories is involved [54] . 
But they are an important ingredient in macroscopic mod¬ 
els and thus should also be checked in mean-field theories. 

It is to be noted that the nine NMP in table [T] taken 
together are fully equivalent to the nine model parame¬ 
ters in the SHF functional ([^, namely Cp, (7^’“, a, Cp, 
and Cp^. There is a one-to-one correspondence between 
the both sets. This allows to consider the NMP also as 
model parameters. This is a more intuitive way to com¬ 
municate the model parameters and it allow direct com¬ 
parison with other mean-field models which are often also 
fully mappable to NMP, see for example the comparison 
of symmetry energies in |56j . 
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isoscalar ground state properties 


equilibrium density: 
equilibrium energy: 






d E 
dpo A 


= 0 


eq 


E 

A 


eq 


side, the charge formfactor is computed as [35] 

Fc{k) = ^ [FqGE,q + Els,gGM]exp ,(17) 

qe{p,n} 

Fg{k) = [d^r expik-rpg{r), 


isoscalar response properties 


incompressibility: 


= Qpl 




effective mass: 

2m* 



2m 


isovector response properties 


E 

dpi A 


eq 


d E 
Otq a 


eq 


symmetry energy: J 

slope of J: L 

TRK sum-rule enh.: ktrk 


1 d^ E 

2 dp\ A 


3 d d^ E 
2 dpQ dp\ A 

2m d E 


h? dri A 


eq 


surface properties (semi-infinite matter) 
surface energy: Osurf see [S3] 

surface symm. en.: asurf.sym see [54] 


Table 1: Definition of the nuclear matter properties 

(NMP). All derivatives are to be taken at the equilibrium 
point of symmetric nuclear matter. 


Pure neutron matter is also an important system as 
there exists an actual realization in neutron stars. Thus 
we will look occasionally at the neutron equation-of-state 
(EoS) E/A {p) which is key input to the computation 

neut 

of neutron star properties, see figure [^ A way to char¬ 
acterize the EoS by one relevant number is to look at the 
slope of the EoS at a typical density. This is d/dpE/A 

neut 

taken at p = 0.1 fm“^. This observable shows up amongst 
others in figure 14 which shows that this observable (as 
most of the neutron EoS) is strongly related to static 
isovector response. 


3.2. Finite nuclei 

The most prominent observables described by an energy- 
density functional are, of course, energy and density. Total 
binding energy is computed in straightforward manner 
with eq. (§. Although the local densities are imme¬ 
diately available as ingredients of the mean-field calcula¬ 
tions, its relation to measurement is more involved and 
deserves some explanation. 

From the experimental side, only the charge density is 
easily accessible through elastic electron scattering which 
allows a more or less model free determination of the nu¬ 
clear charge formfactor Fc(k) [57]. From the theoretical 


where is the form factor of V • Jq augmented by a 
factor /iq/4m^ with pg being the magnetic moment of the 
nucleon, GE,q is the electric form factor and Gm the mag¬ 
netic form factor of the nucleons (assumed to be equal 
for both species). The overall exponential factor takes 
into account the center-of-mass correction for the form¬ 
factor complementing the corresponding energy correction 


(111. It employs the same variance of the total momen¬ 
tum (Pcm) its physical interpretation is an unfolding 
of the spurious vibrations of the nuclear center-of-mass 
in harmonic approximation [40] . The nucleon form fac¬ 
tors GE,q and Gm are taken from nucleon scattering data 

[551 ISO] , for details see m- 

Exact DFT should, in principle, provide a reliable de¬ 
scription of density distributions, or formfactors respec¬ 
tively, in all details [22]. However, actual functionals em¬ 
ploy analytically simple forms which are smooth functions 
of the densities as motivated by a local-density approxi¬ 
mation. It has been shown that this limits the predictive 
value to the regime k < 2k-p in the formfactor where fcp is 
the Fermi momentum |61j . Fortunately, the most crucial 
bulk properties of the nuclear density profile are deter¬ 
mined at low k. These are m 


charge r.m.s. radius: 

3 d^ 

rc = 


Fc{k) 


Fc(0) 

charge diffraction radius: 


fc=0 


Rc = 


4.493 

Tor 


Fc(A:W) = 0 


(18) 


charge surface thickness: 


(TC = 


where kj^'^ is the first zero of the formfactor Eq. The 
diffraction radius Rc parametrizes the overall diffraction 
pattern which resembles those of a filled sphere of radius 
Rc m- The actual nuclear formfactor decreases faster 
than the box formfactor Fbox due to the finite surface 
thickness ac of nuclei which is thus determined by compar¬ 
ing the height of the first maximum of the box equivalent 
formfactor and of the mean-field result Fc- 


log 


Fbox(fc) = 3 


Fbox(^m) 

Fc{km) 
jijkRc) 
kRc 


h — 


5.6 


Binding energy Eb and the three bulk parameters (18) 


of the charge formfactor are the key observables in the 
pool of fit data, see section [4(1] Besides these bulk prop¬ 
erties, we also include in the fit pool odd-even staggering 
of binding energies to fix the pairing strength. It reads for 
neutrons 

AI/HZ,N) = ^iEB{Z,N+l)-2EBiZ,N) + EB{Z,N-l)) 

(19) 
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and similarly for protons. It is interesting to note that 
the odd-even staggering also refers to binding energies. 
But taking the difference of them filters a much different 
information than the one seen in binding energies as such. 
Moreover we consider for the fits a few selected spin-orbit 
splittings in the spectrum of s.p. energies Asni =s i — 

nLj cy 


e ^ 1 to fix the strength of the spin-orbit term. 

Binding energy is also considered for nuclei out¬ 
side the fit pool as predicted (or extrapolated) observables, 
e.g., when estimating the properties of exotic nuclei and 
super-heavy elements. Very useful are also differences of 
binding energies as two-nucleon separation energies or the 
energy expense for a-decay 


Q^ = Eb{Z,N)-Eb{N-2,Z-2)-Eb{2,2) . (20) 


Differences act as amplifying glass. They subtract smooth 
bulk properties and can reveal other aspects of a model as 
mentioned already in connection with the N). We 

will see that also when discussing statistical correlations 
between observables in section 1631 

The charge distribution is mostly sensitive to the proton 
distribution. Unfortunately, experimental information on 
neutron or mass density is plagued with model dependence 
[62] . Promising measurements of neutron radii with high- 
energy particle scattering are coming up [531IM] ; but have 
yet to become more precise. It is worth the effort because 
reliable information on neutron radii will be of invaluable 
help to assess isovector properties of the nuclear functional 
[65] . Thus we will have a look also at the neutron radius 
Tn in the following. We do this in terms of the neutron 
skin from r.m.s. radii 


^skin — ■ (^ 1 ) 

Taking the difference of two similar bulk properties acts, 
again, as an amplifying glass which filters particularly the 
isovector properties. 

SHF is also adapted to deal with nuclear excitations. 
The most prominent ones are the much celebrated giant 
resonances. They are usually described by the random- 
phase approximation (RPA), a dynamical self-consistent 
mean-field theory of small-amplitude oscillations see, e.g., 
[33] . We treat it here with the efficient operator tech¬ 
niques of [551157]. RPA delivers the full spectral strength 
distribution. For the analysis in this paper, it is prefer¬ 
able to have one number to characterize an excitation. 
Giant resonances in heavy nuclei are the perfect candi¬ 
dates for that because their spectrum is well concentrated 
in one resonance. Thus we will consider in the following as 
a measure for typical resonance excitation properties the 
peak energies of the isoscalar giant monopole resonance 
(GMR), the isovector giant dipole resonance (GDR), and 
the isoscalar giant quadrupole resonance (GQR), all in 
208pb q'jje GMR and GQR are very well concentrated 
and thus can be safely computed by the inexpensive fluid- 
dynamical approximation m- The GDR has a somewhat 
broader, though still peaked, distribution and requires full 
RPA for correctly placing the resonance peak. We take 
here the point of view that the giant resonances are well 
described by time-dependent mean-field theory at the level 
of RPA. We have to mention, however, that it is still a mat¬ 
ter of debate whether this description suffices. There are 


calculations including complex configurations (coupling to 
phonons) which indicate that such many-body correlations 
modify the giant resonances [53[. Fortunately, they do 
that mainly for for the width of the resonance. But some 
effects on the peak position cannot be excluded. Another 
problem is that the trend of peak positions towards lighter 
nuclei is not correctly reproduced by SHF-RPA [55]. It is 
not yet clear whether this points to a weakness of mod¬ 
eling isovector density dependence or whether it becomes 
another argument in favour of many-body correlations in 
the giant resonances m- Thus there are still open ques¬ 
tion in the theoretical description of giant resonances. At 
present, we adopt the view that RPA provides relevant 
peak positions, at least for heavy nuclei. We will thus 
consider the resonances in ^°®Pb. 

There is another crucial number which can be extracted 
from the dipole strength distribution, the isovector dipole 
polarizability ao- Just recently, there came up very care¬ 
ful experimental evaluations which deliver ao with high 
precision [ZB [Tg. This is thus a welcome data point 
to assess static isovector response in mean-field models. 
The polarizability can be evaluated as integral over the 
dipole strength with inverse energy weight. This is the 
standard way in experimental evaluation and an option 
in theoretical calculation. Mean-field theory allows an 
even more robust access as the static response to an ex¬ 
ternal dipole field. To that end, one computes the nuclear 
ground state with a small additional dipole field I4xt = 
where d is the isovector dipole operator. The polarizabil¬ 
ity is then deduced from the dipole response of the system 
OiD = d/d\ (d). In our calculation, we use this robust and 
inexpensive option. 

We will also consider fission barriers in super-heavy el¬ 
ements. The fission path is computed by quadrupole con¬ 
strained SHF as a series of mean-field states ‘I*ct 2 o with sys¬ 
tematically increasing quadrupole momentum a 2 o- This 
yields the collective potential for fission as V(a 2 o) = 
A^SHF(‘f’a 2 o)- If to include the full correlation energy 
( [To| . The fission barrier is then 

Bf = Vmax — Vmin + ^'coll.O (22) 

where Vmax is the potential at the peak of the barrier, Vmin 
at the minimum, and Acoii.o is the energy of the collective 
ground state above Vmin- The i?con,o defines the entry 
point for fission as it is explored in experiment. Note that 
we thus include for this observable some correlation effects 
beyond SHF which means that 5/ is not a pure mean-field 
observable. For details of this rather involved definition 
and calculations see [73]|71[. An extensive discussion of 
fission properties in connection with statistical analysis is 
found in m- Here we spot only one example. 

We will also have a look at low-lying quadrupole vibra¬ 
tions (lowest 2+ states). These are related to large am¬ 
plitude motion in mid-shell nuclei and thus treated very 
similar to fission. One generates a collective path with 
systematically changed quadrupole deformation thus map¬ 
ping the collective potential V(q; 2 o) for quadrupole motion. 
The potential is augmented by a collective mass to be com¬ 
puted by self-consistent cranking mass (often called AT- 
DHF mass) and correction for the spurious zero-point en¬ 
ergy. This together defines a collective Schrodinger equa¬ 
tion whose solution then yields the collective excitation 
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energies and associated B(E2) values. For details of this 
rather involved formal framework see, e.g., mm- 


4. Calibration of the model and analysis of predic¬ 
tive power 

4.1. Least-squares fit 

We have used in section [2.3| some formal reasoning to mo¬ 
tivate the form of the SHF functional (§. The actual 
values of the parameters remain open. Up to now, 

there exists no quantitatively successful derivation of the 
parameters from ab-initio calculations and this will remain 
so for a while because the case of nuclear many-body the¬ 
ory is not yet unambiguously settled. It is standard prac¬ 
tice since decades to calibrate the SHF functional (and 
competing mean-field models) to empirical data. Various 
strategies have been used for that in the past. The most 
systematic approach is probably a least-squares (y^) fit 
m which was used in an SHF context first in m and 
has meanwhile become the standard method for calibra¬ 
tion of self-consistent mean-field models. There remains 
a great variety in the bias set by the choice of empiri¬ 
cal data for a fit. To name a few examples out of many: 
some concentrate on spherical nuclei with negligible corre¬ 
lation effects [35], others are particularly concerned with 
deformed nuclei izaiTiiiHni, still others try to adjust also 
spectra of s.p. energies [8T1|82|. The properties of the re¬ 
sulting functionals have much in common, although some 
differences may be found in details. This indicates some 
robustness of the modeling. We will discuss in section [6^ 
the impact of variations of fit data. The standard pool of 
fit data for this publication will be explained at the end 
of this section. Before that, we briefly summarize the 
technique. 

Center piece of y^ fits is a global quality measure by 
summing the squared deviations from the data as 

E 2 

Xobs 5 

obs 

^(th) _ ^(exp) 

E '^obs,nucl '^obs,nucl /OQ^ 

Aabs.nucl ’ ^ ^ 

nuci 

where O stands for an observable, “obs” for a type of 
observables {E, r^^s, •■•)> “nucl” for a nucleus (defined 
by Z,N), the upper index “th” for a calculated value, 
and “exp” for the experimental value. The denominator 
AOobs.nuci quantifies the adopted error of that observable. 
It renders each contribution dimensionless and regulates 
the relative weights of the various terms. It is to be noted 
that AOobs.nuci cannot be identified with the experimen¬ 
tal error on the given observable which is usually much 
smaller. The limiting factor are limitations at the side of 
the model, i.e., the quality we can expect from a mean- 
field description. A self-regulating choice is to tune all 
AOobs.nuci in one group “obs” the same way such that 
Xobs yields in the average a contribution of about one for 
each nucleus [76l |83| , for details in the nuclear context see 
[SI]. 

The total quality measure is a function of all model 
parameters, i.e. y^ = y^(p) where p = {pi, The 


X^ = 


Xobs 


optimal parameters Pq are those which minimize y^, i.e. 

Po ■ X^iPo) = Xo = min ■ (24) 

The minimium condition looks straightforward. But find¬ 
ing the absolute minimum is a non-trivial task because the 
y^ landscape has several minima and often discontinuities 
in between |26j . We employ the minimization technique 
from Bevington m which works well in the vicinity of 
minima and combine that with Monte-Carlo steps to ex¬ 
plore the rongh y^ landscape in a broader range. For large 
scale searches, there are more elaborate methods around, 
e.g. more robust iteration to a minimum [85] or genetic 
algorithms for particularly obstinate cases [86] . 

Finally, we comment briefly on the choice for the pool of 
fit data used later on. Some more details are provided in 
appendix]^ We are using exactly the same data as in the 
survey of [32]. In that paper, the experimental values and 
the adopted errors are explained in great detail. Thus we 
need not to repeat it here. The set of fit data includes the 
bulk properties rmis, .Rdiffr, and a plus pairing gaps 


and some spin-orbit splittings (see section 3.2). The nuclei 


are restricted to be spherical and carefully selected to have 
small ground-state correlations [41] . This guarantees that 
we deal with nuclei which can reliably well be described 
by a mean-field model. Another group which has small 
correlation effects would be strongly deformed nuclei. We 
omit them for practical reasons because they are costly to 
compute. All transitional nuclei require corrections from 
vibrational correlations and must be discarded from any 
fit pool. 


4.2. Statistical error analysis 


Not only the optimal parametrization pq, but also the pa¬ 
rameters in the vicinity of po deliver a reasonable repro¬ 
duction of data. This is systematically quantified in sta¬ 
tistical analysis [TH] [83] . Although we will see later that 
the residual errors 0^b?„uci " ^’^ud (also called sim¬ 
ply “residuals”) in y^ fits of nuclear mean-field models are 
not really statistically distributed, we will employ the well 
developed schemes of statistical analysis for estimating ex¬ 
trapolation errors and correlations between observables. It 
remains in any case a powerful tool to explore the y^ land¬ 
scape and the information it contains about the interplay 
of model and fit data. 

Assuming a statistical distribntion of residuals, one pos¬ 
tulates a probability distribution of reasonable model pa¬ 
rameters as [13 [ST] 

IU(p) oc exp(-y2(p)) . (25) 

Their domain is characterized by y^(p) < yo + 1 (see 
Sec. 9.8 of Ref. [13)- The range of p fulfilling this condi- 
tion is usually small and we can perform a Taylor expan¬ 
sion 


X^{p) 


Np 

Xo + X! iPa-PO,a){C~^)a0(,P0-PO,0),{‘2‘Q) 

a,/3—1 


— 3 3 

2 ^Pa ^p& A- 




Jif-. 


(27) 


dp.O, 


AO, 


( 28 ) 







where J is the rescaled Jacobian matrix and C the covari¬ 
ance matrix. The latter plays the key role in covariance 
analysis. The domain of reasonable parameters is thus 
given by 


X {p)-Xo~P-C • p < 1 


(29) 


which defines a confidence ellipsoid in the space of model 
parameters. It is related to the Taylor expanded proba¬ 


bility distribution (251 which becomes 


W{p) = (27r)^det{C} 


- 1/2 1 
exp(--p 


P) (30) 


where the fore-factor guarantees proper normalization 

^d^pW = 1. 

Any observable A is a function of model parameters 
A = A{p). The value of A thus varies within the confi¬ 
dence ellipsoid, and this results in some uncertainty AA. 
Usually, one can assume that A varies weakly such that 
one can linearize it 


A{p) 


A{po) + G^ -{p- po) 
dpA 


(31a) 

(31b) 


This together with the Gaussian probability distribution 
(30) allows to compute analytically averages and vari¬ 
ances. The average becomes 


A = J d^pW{P) A{p) = A{po) (32a) 


5. Strategies for estimating errors 

Although very powerful, a fit is a black box and se¬ 
duces the user to use it without much understanding of the 
physics beyond. One plugs in a model, chooses a couple 
of relevant fit data, and grinds the mill until one is con¬ 
vinced to have found the absolute minimum Xo together 
with the optimal parameters Pq. Thus far one may find 
a satisfying description of the fit data but may live with 
little understanding of the relevance and reliability of the 
model. The challenging task remains to understand the 
model thus achieved, in particular its reliability in extrap¬ 
olations to other observables. This is the quest for error 
estimates, see [H] for a basic discussion of error estimates 
in the context of nuclear models. The problem does not 
have a simple and unique answer. Statistical errors as 
explained in section |4.2| are well under control. The hard 
part is to estimate the systematic errors coming from prej¬ 
udices in the modeling or choice of data. The bad news 
is that there is no systematic way to assess systematic 
errors. From a critical point of view, every physical the¬ 
ory holds only preliminarily and is always threatened by 
hidden insufhcencies which may be revealed some days by 
new observations. As a consequence, we need to handle 
our theories with a great deal of restless error awareness. 
The best way to enhance confidence in a model is to scru¬ 
tinize it employing a great variety of different checks and 
so piece-wise put together a view of the sources of uncer¬ 
tainty. In the following, we summarize a couple of possible 
checks: 


as one would have expected. The variance quantifies the 
fluctuation of A around A{po) and is 


AA = \/a^ , (32b) 

= j d^pW{P),{A{p)-Af 

OL^ 


with as given in eq. (31bI. As it is derived from 
a statistical interpretation of the landscape, it is also 
coined statistical error. Along the same lines, one can also 
define a cross-variance between two different observables. 
It reads 


AAAB 


J d^pW{P), (A(p) - A){B{p) - B) 
Y,GtCo.pG$ (32c) 

OL^ 


This then allows to defined the covariance, or correlation 
coefficient, 


I AAABI 
AAAB 


(32d) 


It quantifies the statistical correlations between two ob¬ 
servables A and B. A value cab = 1 means fully cor¬ 
related where knowledge of A{p) fully determines B{p). 
A value cab = 0 means uncorrelated, i.e. A{p) and B{p) 
are statistically independent where knowledge of A(p) car¬ 
ries no information whatsoever on B{p), for examples see 

[SHI [89]. 


1. Exploring a model from the theoretical side 


1. Estimates of beyond mean-field effects: 

The first step at all is to check whether the 
model under consideration is appropriate to de¬ 
scribe the wanted data. This means for SHF 
to check many-body correlations beyond mean 
field. Density matrix expansion (see [121 150] 
and sectior2.3) indicates that short-range cor¬ 
relations can be mapped into an energy func¬ 
tional. However, correlations from low lying 2'^ 
states fluctuate strongly and cannot be mapped 
into a smooth functional. Nuclei which acquire 
large correlation corrections should be excluded 
from the pool of fit data [32l [4T| . The search for 
possible correlation effects has not yet come to 
an end, for an example see section [6.4| 


2. Exploring the stability of the model: 

Models which rely on expansions are likely 
to have regions of instability. There is also 
the problem that x^{p) can have non-analyical 
points. It is crucial to explore those dangerous 
regimes and set rules to avoid them, or to cure 
the problem. These considerations have accom¬ 
panied the development of SHF all along, for 
recent critical studies, particularly of the tensor 
term, see pSl[?7l 1^[^ . 


3. Formal analysis of inter-dependencies: 

Correlation analysis (see point 2.2) can reveal 
dependencies between observables. To make 
sure that these are more than statistical corre¬ 
lations, one ideally establishes the connection at 
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a purely formal level. Literature on such stud¬ 
ies is unoverseeable. We mention here only two 
examples, one for the relation between effective 
mass m*/m and level density [90] and another 
one for the mapping of SHF parameters to NMP 
via the LDM pT] . 

4. Variations of the functional: 

A test for a model from within is to vary the ba¬ 
sic ansatz. This means, e.g., for the SHF model 
to omit or to add terms or to vary the den¬ 
sity dependence [SI]. Another way is to com¬ 
pare and/or accumulate the results from differ¬ 
ent models as done by combination with RMF in 


tical error on extrapolations. It serves as useful 
indicator for safe and unsafe regions of the model 
see, e.g., section 

2. Correlations between observables: 


The covariances (32d) reveal dependencies be¬ 


groups of observables “obs” as, e.g., energy or 
radius. One can omit a group from the pool of 
fit data. Comparison with the full fit allows to 
explore the impact of the omitted group, see [M] 
and section Hill 


4. Test of predicted observables: 

A natural test is to compare a prediction (ex¬ 
trapolation) with experimental data wherever 
available. It is a strong indicator for a systematic 
error if the deviation exceeds significantly the ex¬ 


trapolation error (32b). An example for such a 


|89l|93]. We will exemplify in section 6.1 adding 
and omitting terms. 

2. Straightforward statistical analysis 

1. Extrapolation uncertainties: 

This employs eq. (32b) to estimate the statis¬ 


critical case is the energy of super-heavy nuclei 
where systematic deviations point to a problem 
of SHF, see [53] and section |6.0.11 


tween observables. This helps to estimate the 
information content of a new measurement, for 
examples see [55] and section 6.5 It may be 


combined with a variation of fit data (see point 


3.3) to distinguish between model properties and 


data dependence. It can also serve to motivate 


searches for formal connections (see point 1.3). 

3. Sensitivity of model parameters p: 

The Jacobian matrix J together with the covari¬ 
ance matrix C allows to explore the impact of 
each single model parameter pa on a given fit ob¬ 
servable Oi- Examples are found in 

3. Exploring the model by dedicated variations 

1. Unresolved trends of residual errors: 

A perfect model should produce a purely 
statistical (Gaussian) distribution of residuals 

C’Snuci - ^Suci- Well visible trends of resid¬ 
uals of one sort of observable over the series of 
nuclei indicate deficiencies of the model. All nu¬ 
clear mean-field models produce still strong un¬ 
resolved trends see, e.g., [32E11|SS| and section 

Eon 

2. Dedicated variations of parameter or observable: 
As a complement to correlation analysis (see 
point 2.2), one obtains a direct view of depen¬ 


The above list of testing strategies shows that establishing 
a model requires a broad palette of tools and counter¬ 
checks. The present paper can not exemplify all of them. 
We will address the groups and [^ as well as point 1^ 
from groupie 

6. Results 

6.0.1. Bulk properties 

Modern SHF parametrization have reached a high quality 
of reproduction of nuclear ground state properties. As ex¬ 
ample, we quote here the r.m.s. errors for SV-min which 
was fitted to ground state properties of semi-magic spher¬ 
ical nuclei (32). They are 0.6 MeV for binding energies 
Eb, 0.029 fm for charge diffraction radii 0.022 fm for 

charge surface thicknesses, and 0.017 fm for charge r.m.s. 
radii. Similar quality is found for all modern parametriza- 
tions which appeared in the last decade. 

Although small, the question is how the errors are dis¬ 
tributed over the nuclei. We do that here for SV-min [32] . 
Figure]^ shows the residuals of ^.m.s. 

radii (upper panel) and of binding energies (middle and 
lower panel) for all nuclei for which data were available. 
Different symbols (and colors) distinguish nuclei from the 
fit pool, well deformed nuclei, and everything else. The 
lower panel shows binding energies Eb from SHF calcula¬ 
tions with c.m. correction, but without any other correla¬ 


tion correction (10). The fit nuclei gather nicely within a 


dencies by dedicated variation of a parameter 
Pa or observable A{jp). To that end, one fits a 
parametrization with constraint on a fixed value 
of Pa or A and systematically varies the con¬ 
strained value [32] ■ Then one looks at the trend 
of other observables B{pa), or B{A) respectively. 
This test is particularly instructive when using 
variations of NMP, see [32] and section 6.3 


3. Variations of fit data: 

The fit observables Oobs nuci stem from different 


small error band as given by the above mentioned r.m.s. 
errors. That is no surprise because they were selected to 
have negligible correlations [321 El] • The other nuclei devi¬ 
ate significantly. The middle panel shows the results with 
all corrections from collective correlations, see eq. ( [l0| . 
This brings all results closer to a narrow error band. The 
deformed nuclei in the medium heavy region {A < 200) 
perform already very well while transitional nuclei still 
overshoot. 

The errors are obviously not statistically distributed 
but show strong trends which are not resolved by the 
given SHF model. Even with correlation corrections (mid¬ 
dle panel) errors are smallest towards shell closures and 
largest in between, particularly for the transitional nu¬ 
clei (red triangles). This is clearly a shell fluctuation not 
fully accounted for by the present estimate of correlations 
(only from the lowest 2+ state). More than that, there 
is a strong trend to underbinding with increasing system 
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Fig. 1: Deviations of SHF calculations from experimental data 
for all nuclei for which data were available. Fit nuclei are shown 
with filled blue squares, well deformed nuclei with green circles, 
and all other with red triangles. Isotopic chains are connected 
by a line. Lower panel: for binding energies computed with 
(deformed) SHF including only the c.m. correction. Middle 
panel: for binding energies with SHF plus corrections from 
collective correlations (vibration and rotation). Upper panel: 
For charge r.m.s. radius computed with (deformed) SHF. 


size A for deformed nuclei (green circles). This becomes 
obvious in the regime of actinides and super-heavy ele¬ 
ments (A>220). Even attempts with variations of density 
dependence |92j and including deformed nuclei in the fit 
[55] could not resolve the problem. It seems to be in¬ 
herent in the way the present SHF functional ([^ models 
density dependence. For example, traditional RMF func¬ 
tionals which model density dependence quite differently 
by non-linear self-coupling of the tr meson have the op¬ 
posite problem: they tend to overbind super-heavy nuclei 
[55] . The unwanted trends in figure suggest that there 
is still a need for improvement in nuclear self-consistent 
mean-field models. 

The upper panel shows deviations for charge r.m.s. 
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Fig. 2: Evolution of errors for a chain of super-heavy elements 
with increasing proton number Z in steps of 2. The neutron 
number has been increased accordingly in steps of 4, occasion¬ 
ally steps of 2, to cover the landscape of super-heavy elements 
ending up near the magic neutron number N = 184. Shown 
are extrapolations errors (denoted “extrap.”) and deviations 
from data Esuf — Aexp (denoted “deviat.”) where Ashf in¬ 
cludes correlation corrections as explained in connection with 
figure [^ The parametrization SV-min was used. 


radii. Correlation effects are smaller for this observable 
and not shown. There are also less data which leads to less 
points on the plot. The deviations gather better within 
the error band from the fit nuclei. Nonetheless, they are 
also not statistically distributed. We observe again the 
shell fluctuations between shell closure and mid-shell nu¬ 
clei. The pronounced trend for deformed, super-heavy nu¬ 
clei which was seen for binding energies is absent for r.m.s. 
radii. 


6.0.2. Far extrapolations 


An important quantity delivered by the technique is 
the variance (32b) of an observable. It provides an esti¬ 


mate for the extrapolation uncertainty which is extremely 
useful to get an impression how reliable an extrapolation 
may be. We exemplify that here for an extrapolation to 
super-heavy nuclei and an extremely far extrapolation to 
neutron star properties. The case of super-heavy nuclei 
continues the discussion from the previous section |6.0.I| of 
an unresolved trend for heavy nuclei. 

Figure [2] sho ws extrapolation uncertainties computed 
with eq. ( |32b| ) and actual deviations from data (where 
available) computed with SV-min for a chain of super¬ 
heavy nuclei from Z = 90 to Z = 120. The sequence of 
neutron numbers has been chosen to cut a path through 
the landscape spanning from Z/N=90/126 to 120/182. It 
was: N =126, 130, 134, 138, 142, 146, 148, 150, 152, 154, 
158, 162, 166, 170, 174, 178, and 182 The extrapolation 
uncertainties (shown as error bars) grow steadily when 
moving away from the region of fit nuclei. This is an ex¬ 
pected behavior: the farther away the extrapolation the 
larger the uncertainty. The chosen sequence still allows 
to compare with data for more than half of the sample. 
The actual deviations from available experimental values 
(shown as filled circles) are significantly larger than the 
estimates of uncertainty from statistical analysis. The fig- 
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Fig. 3: Mass-radius relation of a neutron star computed for SV- 
min (heavy blue line). The shaded area indicates the range of 
the extrapolation errors. The mass-radius correlation ellipsoid 
is shown by a read line. The dashed horizontal line indicates 
the mass of the heaviest identified neutron star. The black dots 
with errors stand for the maximum mass and maximum radius 
from SV-min. Data taken from m- 


ure demonstrates nicely the possible devastating effect of 
a systematic error. In this case we are in the lucky situ¬ 
ation that we could already identify the problem because 
there are sufficient data for checking, see figure and the 
discussion thereof. The case becomes uncontrollable for 
totally new observables for which one does not dispose of 
benchmark data. Nonetheless, with the statistical extrap¬ 
olation error we have an indispensable tool to get an idea 
about the uncertainties to expect and a safe lower esti¬ 
mate for them. Mind the the extrapolation uncertainty in 
figure has at least the correct order of magnitude. 

Figure addresses as next example the extrapolation 
to neutron star mass and radii. The heavy line shows 
the results of the Tolman-Oppenheimer-Volkov equation 
for the mass-radius relation using as input the neutron 
equation-of-state from SV-min [96) . The faint horizontal 
line indicates the maximum mass of a neutron star found 
so far m- This sets a lower limit for the maximal mass 
because one cannot exclude that an even heavier star may 
be found some day. The prediction from SV-min stops far 
below the wanted maximum mass. However, the outlook 
improves when taking into account the extrapolation er¬ 
rors, marked by the shaded area. The uncertainties are 
huge which is understandable because we deal here with 
a really far extrapolation, from finite nuclei to compar¬ 
atively infinite star matter. The good news is that the 
error band scratches the wanted maximum mass. This in¬ 
dicates that inclusion of wanted neutron star properties 
can be accommodated by a fit without sacrifices. This 
was done with success in |96) yielding the parametrization 
TOV which indeed hits correctly the mass line and has 
also a much smaller uncertainty band due to the inclusion 
of star data. The figure contains one more detail. This 
is the error ellipsoid plotted inside the error band. It indi¬ 
cates the region in which the probability distribution (30) 
is larger than 1/2, this is the region of predictions from 
the range of “reasonable parametrizations”. In this exam¬ 
ple, the ellipsoid looks more like a needle going diagonal 
through the error band. This means that we encounter 


here a highly correlated scenario where a given value of 
radius fixes inevitably a corresponding value for mass. It 
is a prototype example for a case with large covariance 
(|32dl close to 1. 


6.1. Hierarchy of terms in the functional 


The SHF functional (|^ was motivated by a Taylor expan¬ 
sion of the a supposed microscopic effective interaction, 
see section |2.3| This implies a hierarchy of importance of 
the terms in the functional with the zeroth order of ex¬ 
pansion (the purely density dependent parts) representing 
the leading terms. A further ordering principle comes from 
the fit data. The landscape of known nuclei extends widely 
over size A but only over a narrow band of isotopic chains. 
This means that isoscalar properties are well determined 
by the long trends with A while knowledge of isovector 
properties is limited by the small extension in isovector di¬ 
rection N — Z together with the fact that only one sort of 
density is well known empirically, namely the proton den¬ 
sity measured by electron scattering I^. This combined 
sorting principles led to the assignment of a “minimal set 
of terms” as indicated by yellow (gray) shading in eq. ®. 
In this section, we are going to explore the importance of 
terms in detail by starting from the most sparse model 
using only density dependent terms and adding step by 
step more terms up to the full model At each stage, 
the model consisting out of the chosen terms is optimized 
by minimization of always using the same data pool 
from [32], see end of section 4.1 and apendixThis tries 
to make the best out of the given functional. The strat¬ 
egy differs from a similar investigation of stages in [39] . 
There, we started from ab-initio calculations in bulk mat¬ 
ter to define the background of LDA similar as in electron 
systems |22j and added stepwise derivative terms with in¬ 
formation on finite nuclei. Here we optimize each stage as 
far as possible. 

We consider the following stages of the model: 


LDA: 

This first stage invokes the terms with the parameters 
Cg, Cf, and a. It covers only the purely den¬ 

sity dependent terms which are active in homogeneous 
matter. We call it LDA because that choice produces 
a purely density dependent functional (no derivative 
terms wherever). But we optimize it empirically. It 
is not LDA in its strict sense where exact solution for 
bulk matter are used to deduce a functional 


-I-grad: 

To the previous stage, the term oc is added. This 
gives the model the flexibility to account for surface 
properties. 


-I-Is: 

To the previous stage, the isoscalar spin-orbit 
term is added. This allows to put shell closures at 
the right place. 

+Cf’“: 

To the previous stage, the density dependent isovec¬ 
tor term oc (7/’“ is added. This allows to resolve 
more isovector trends in the ground state properties, 
if there are still any. 
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Fig. 4: Average and r.m.s. errors on basic bulk properties, 
binding energy Eb, charge r.m.s. radius rrms, charge diffraction 
radius iZdiffr, and charge surface thickness Uaurf- Averages are 
taken over the pool of fit nuclei. The average deviation from 
data is shown by filled boxes and the r.m.s. deviation by the 
error bars around the boxes. Results are shown for the series of 
fits with increasing number of terms in the functional § , see 
text. The lower panel for errors on binding energy is split into 
two energy scales. The first two entries (“LDA” and “+grad”) 
are plotted with respect to the scale on the left y-axis, the 
other three with respect to the right y-axis which is stretched 
by factor 5 with respect to the left scale. 


SV-min: 

This is result from invoking the full SHF functional 
Q. It was published in [32] an discussed in great 
detail there. 

In all cases, we include, of course, the pairing functional 
Q. The first two stages (“LDA” and “-l-grad”) are too 
rough to deserve a detailed pairing adjustment. Here we 
use volume pairing, i.e. po.pair —t oo, with fixed Ipair,p = 
Ipair.n = 300 MeVfm^. For the further three stages which 
provide a detailed modeling, we optimize all three pairing 
parameters together with the free parameters of the SHF 
model. 

Figure l^shows for the above listed sequence of stages of 
the model average and r.m.s. deviations for the four fitted 
bulk properties where averages are taken over the pool 
of fit data. The errors for the simplest stage, LDA, are 
large, particularly for the binding energy where the error 
band reaches ±20 MeV. But this is, in fact, a satisfying 
result. It shows that already these few leading terms are 
capable of accommodating roughly the energies over all 
A, and even better for radii and surface thickness. Of 
course, details are bound to be wrong because we are yet 
missing the spin-orbit term and thus put shell closures at 
the wrong nucleon numbers. 

Adding the gradient term allows to tune the nuclear sur¬ 
face energy. This improves the average quality of energy 
by half an order of magnitude. With ±4 MeV it is as far 
as on can come without appropriate shell effects. The de¬ 
scription is even better for radii and surface thickness be¬ 
cause shell fluctuations play a minor role in these observ¬ 
ables. Here one has reached almost the final quality. The 
great improvement due to the gradient term shows that 
the surface energy is a crucial contribution to a nuclear 
energy functional. This complies with the same findings 
in the purely macroscopic droplet model of the nucleus 

mm- 

The next substantial improvement is achieved with 
adding the spin-orbit force, stage “±ls”. This yields a 
mean-field model where shell closures emerge at the right 
nucleon number and also all other shell effects acquire the 
correct place and magnitude. Consequently, we see an¬ 
other dramatic reduction of the energy error. The model 
has now reached a high quality. All further terms deliver 
only gradual improvement. The stage “±ls” with its 6 free 
SHF parameters thus represents a minimal model which 
produces already an excellent description of the fit observ¬ 
ables (bulk properties of the nuclear ground states). 

In the step to stage “±C'f’“”, we explore the gain by 
adding more isovector flexibility. There is minimal im¬ 
provement in the quality of radii and surface thickness and 
there is a 25% reduction in the energy error, not dramatic 
but non-negligible either. There is thus sufficient isovec¬ 
tor information in the fit pool to fix this second isovector 
parameter. 

Finally with “SV-min”, we have activated all terms in 
the functional ([^. The gain as compared to the previ¬ 
ous stage is extremely small. The additional terms in 
that step are all of isovector type and the pool of fit data 
does not contain enough isovector information to fix them. 
Thus “SV-min” leaves large uncertainties in these isovec¬ 
tor parameters. Further data are required to determine 
the isovector model parameters more precisely. Isovec- 
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Fig. 5: Trends of deviations from data for the series of fits 
with increasing number of terms in the functional, see section 

The lower panel shows deviations for the binding energy Eb 
and the upper panel for the charge r.m.s. radius rrms- Results 
along isotopic or isotonic chains are connected by a line. 

tor response properties (GDR, dipole polarizability, and 
NMP, see section [63| ) and the neutron skin (21) add valu¬ 
able information as we will see later. 

Figurej^shows the detailed map of deviations from data 
for the above introduced stages of the SHF functional. For 
simplicity, the plot is limited to errors in energy and r.m.s. 
radius. The diffraction radius and surface thickness be¬ 
have similar. The errors for LDA show a strong trend over 
mass number A. The fit for LDA has found a global com¬ 
promise between light and heavy nuclei. But the model 
at LDA level is not capable to deal properly with system 
size. The step to “-l-grad” resolves the global trend at 
once. A straight line through the errors is identical with 
the zero line. But there remain the strong shell fluctu¬ 
ations which are due to the fact that shell closure lie at 
the wrong nucleon numbers. These fluctuations are much 
reduced with the step to “-|-ls” whose energy errors (lower 
panel) then gather closely around the zero line. The errors 
on radii (upper panel) seem to be more resistive. But this 
impression is to a large extend a matter of scale. LDA 
delivered already quite good values for the radii. In turn, 
the discrepancy between LDA and further stages is not as 
dramatic for the radii as it is for energy (see also figure |^. 
Thus we see the deviations on radii with higher resolution 
and that reveals more clearly the irresistible unresolved 
trends. 

The further steps deliver only gradual improvement, 
hardly distinguishable at the scales of this plot. The 
remaining errors, however, are still not statistically dis- 
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Fig. 6: Predictions for giant resonances and dipole polariz¬ 
ability in ^°®Pb compared with experimental data. Results are 
shown for the series of fits with increasing number of terms 
in the functional. A new entry as compared to figure and 
l^is the parametrization SV-bas. It is fitted the same way as 
SV-min but includes constraint an NMP in order to allow also 
a good reproduction of giant resonance energies, see m- 


tributed. There remain some unresolved trends seemingly 
related to shell structure (shell closures versus mid shell). 
This calls for further refinement of the model. But it is 
yet open whether this can be achieved within the SHF 
and pairing functional or whether we are forced to invoke 
correlation effects. 

Figure shows the peak energies for giant resonances 
and the dipole polarizability for the stages of the SHF 
functional, all for the nucleus ^°®Pb. The lower panel 
shows the GMR. It improves with each stage and reaches 
already at the minimal model, stage “-|-ls”, good agree¬ 
ment with the experimental value. This is not surpris¬ 
ing. The GMR is closely tied to the incompressibility 
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K which characterizes the response to changes in density 
[981199] . The incompressibility, in turn, depends on the 
parametrization of density dependence in the model. And 
this seems to be well determined by the fit data which 
span a broad range of mass numbers (isoscalar informa¬ 
tion) . We learn that the isoscalar aspects are already well 
converged with the stage “-|-ls”. 

The second panel from below shows results for the GDR. 
They are all far below the experimental point even. The 
reason is clear for all stages before SV-min because these 
do no invoke the elfective mass terms and it is known 
that the GDR depends sensitively on TRK sum-rule en¬ 
hancement ktrk which is closely related to Cl [1001 [sg. 
But even the step to for SV-min, which includes all terms 
of the functional does not improve the situation. The fit 
data are not strict enough to push ktrk to the right value. 
On the other hand, they do not fix ktrk so much. The un¬ 
certainty in the GDR energy is 1.4 MeV for SV-min. This 
leaves sufficient leeway to tune ktrk- This is done with 
the parametrization SV-bas which optimizes some NMP 
to tune the giant resonances in ^°®Pb|32]. This leads ob¬ 
viously to a good reproduction of the GDR. 

A similar situation is encountered for the GQR (second 
panel from above). All stages up to including SV-min 
underestimate its resonance energy. This is no surprise 
for the stages LDA through because they do not 

invoke the term oc Cq and thus have a mean field with 
effective mass m* jm = 1. But the GQR is known to 
depend sensitively on ra* jm = 1 and require a lower value 
0.8-0.9 jlOlj . This option is set free with SV-min. Similar 
as for the GDR, SV-min does not put in*/m, and with it 
the GQR, at the right place, but leaves sufficient leeway 
to tune it properly. This is again achieved by SV-bas. 

The upper panel in figurej^shows the dipole polarizabil¬ 
ity «£) for the various stages of the functional. It starts 
for LDA with a rather wrong value, improves with each 
step, and reaches the data already at the stage 
where the functional was given sufficient isovector freedom 
to take advantage of the (static) isovector information in 
the fit data. 

It is interesting to note that two data points, GMR (re¬ 
lated to incompressibility K) and an, are well described 
already from the fits to ground state data while the two 
other, GDR and GQR, need extra measures. Mind that 
K and ao are static response properties. No surprise then 
that they are well determined by ground state data. The 
GDR and GQR, on the other hand, are related to the 
dynamic response properties to*/to and ktrk- It is plau¬ 
sible that dynamical properties require dynamic data for 
calibration. Thus we need some information from excited 
states to fully fix the functional §. 


6.2. Impact of groups of data 


In this section, we take up point 3.3 of the analyzing strate¬ 
gies outlined in section the variation of fit data. An 
analysis following this strategy was already given in [94j . 
We repeat part of that material for a reduced set of vari¬ 
ations and, on the other hand, add new variations. Basis 
is the pool of fit data as developed in [35], see also end of 
section |4.I| and appendix [^ An unconstrained fit to the 
full set yields the parametrization SV-min. Now we are 
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Fig. 7: Average and r.m.s. errors on basic bulk properties, 
binding energy Eb, charge r.m.s. radius rrms, charge diffraction 
radius Rdiflfr, and charge surface thickness Uaurf. Averages are 
taken over the pool of fit nuclei. The average deviation from 
data is shown by filled boxes and the r.m.s. deviation by the 
error bars around the boxes. Results are shown for the series of 
fits with different selections of groups of data as listed in table 

H 

interested in fits with deliberate omission of groups of fit 
data. This yields the parametrizations as listed in table 



Eb 

f’rms.C 

RdiS,C 

O-C 

“half data” 

x/2 


x/2 

x/2 

SV-min 

X 

X 

X 

X 

“A-brrms only” 

X 

X 

- 

- 

“A only” 

X 

- 

- 

- 

“rrms only” 

- 

X 

- 

- 


Table 2: Included data sets from the standard pool of 
fit data from finite nuclei (33] ■ A “x” means included, 
stands for excluded,and “x/2” means that only half 
of the data were included. Pairing gaps Apair and spin- 
orbit splitting Els are included in all data sets and not 
listed explicitly above. The observables were explained in 
section I3.2I 

We have also studied omission of £is or Apau-. This showed 
only minor effects and will be not reported here. 

Figure [7] demonstrates the effect of varied fit data on 
the average and r.m.s. residuals taken separately over a 
group of data. Mean values deviating significantly from 
zero within the scale set by the uncertainties indicate some 
basic incompatibility of the observable with the model. 
Ghanges on the r.m.s. error indicate the sensitivity to a 
group of observables. 

First, we should realize that the scales in this figure are 
generally much smaller than those of the various stages of 
the functional in figure All fits with reduced data sets 
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produce very similar models. Even when taking “i? only” 
or even “rj-ms only” we obtain an agreeable model. This 
gives us some confidence that the model has predictive 
power and allows extrapolations to other observables or 
other regions of nuclei. 

However, looking closer at the results in figure we see 
at some places differences which carry some information 
on tensions within the model and at other places continued 
similarity which indicates the strongholds of the model. 

Let us first compare the set “half data” with SV-min. 
Both sets produce exactly the same averages. In fact, look¬ 
ing at detailed residuals (not shown here) we find that both 
fits produce very similar results. This can be read off also 
globally from the fact that the r.m.s. errors (taken over 
all data) are similar. This proves that the SHF functional 
interpolates well. 

The step from SV-min to “E+r^yns only” has dramatic 
consequences for i?diffr,c and ac- The average error in 
these two observables makes a large jump upward and 
similarly the r.m.s. error. This indicates that there is 
some incompatibility between and i?diffr,c in the 

present model. That could already be spotted in SV-min. 
Already here, we find a large average error for i?diffr,c- In 
fact, the r.m.s. error is almost exhausted by the average 
error. The same behavior persists in the step to “A only”. 
But here, we also find the same increase in the average 
error for rmis.c which indicates a slight incompatibility of 
Eb and rrms.c- It is likely that both incompatibilities, 
for i?diffr,c as well as for rpms c are related to insufRcen- 
cies in the surface profile and thus to the model of density 
dependence. 

The average error of energy Eb is near zero for all hts 
which include Eb- This indicates that the SHF functional 
is well suited to adjust the energy over all ranges of nuclei. 
The model “rrms only” produces a slight mismatch of the 
average. However, these stay well within the huge error 
bars in this case. The large size of the error bars can be 
viewed as an “automatic” warning from the analysis. 
It signals that the data set “crms only” does not contain 
sufficient information to fix the energies. It is remarkable, 
though, that the other direction, determining r^nis by a 
fit “E only”, works much better. Energies are the leading 
input to the pool of fit data. 

Figure shows the effect of varied fit data on pre¬ 
dicted/extrapolated observables, nuclear matter proper¬ 
ties (NMP) and key observables of the three nuclei ^°®Pb, 
266 Hs and Z=120/N=182. For Z=120/N=182, we con¬ 
sider binding energy Eb and a-decay energy accord¬ 
ing to eq. (20). No data are available yet for this nu¬ 
cleus. For ^‘’*^Hs, we consider binding energy Eb nnu 
and fission barrier Bf [im]. The experimental value for 
Bf is given with error bar as it is associated still with 
large uncertainty. The binding energy includes the rota¬ 
tional zero-point energy (12) which is obligatory for de¬ 
formed nuclei. In ^°®Pb, we consider response properties, 
the peak energies of GMR, GDR,and GQR and the dipole 
polarizability an (see section 3.2). The NMP considered 
are: binding energy E/A, density incompressibility 
K, isoscalar effective mass m*/m, symmetry energy Osym, 
and Thomas-Reiche-Kuhn (TRK) sum rule enhancement 
factor ktrk, all taken at the equilibrium point of sym¬ 
metric matter. Figureshows the effect of variation of fit 


data on the observables and their uncertainties. Changes 
in uncertainty indicate the importance of the omitted data 
group on the observable. A shift of the average shows what 
data are pulling in which direction. 

The effects on NMP (left two columns) are generally 
large. A most pronounced shift is produced by omitting 
the formfactor information of i?diffr,c and ac in “F -|- Crms 
only” and “A only”. This leads to a large jump in bulk 
equilibrium density pnm and incompressibility K. There 
is also a jump in the isovector response Osym in addition 
to the generally strong changes. It is also interesting to 
note that formfactor information keeps the effective mass 
m*/m down to values below 1 while fits without radii let 
m*/m grow visibly above one. The reason is probably 
that m*/m has an impact on the surface profile, thus on 
^'rms.c and ac, and, in turn, also on Rdiffr.c- This sug¬ 
gests that the formfactor is closely related to the surface 
prohle of the nucleus. This, in turn, seems to be con¬ 
nected with modeling the density dependence as one can 
read off from the strong impact on K (second derivative 
with respect to density). The step to “rrms only” restores 
approximately the values of pnm and K from the full ht 
SV-min. This indicates that there is also some competi¬ 
tion between and Eb- These two quantities pull in 
different direction which indicates that they would prefer 
different surface profiles. Mind, however, that the devia¬ 
tions in the average errors are fully compatible within the 
error bars (characterizing the r.m.s. errors). This means 
that all choices of subsets of data are compatible when tak¬ 
ing into account the message about uncertainties of a ht 
which is quite appropriately delivered by statistical analy¬ 
sis. Nonetheless, the observed variations of average errors 
indicate that the SHF functional may deserve further hne 
tuning to reduce these conflicts between observables. The 
most likely aspect to work on is the density dependence of 
the functional. 


The variances usually grow when omitting data which is 
plausible because less data mean less determination. Par¬ 
ticularly large errors emerge for the sets “E only” and 
“j’rms only”. This shows that combined information with 
the charge formfactor is extremely useful to conhne the 
parametrization. But mind that the error bars for “E 
only”, although being generally larger, are still in accept¬ 
able ranges. This shows that data from energy alone can 
already provide a reasonable parametrization. One may 
argue similarly with the errors for ‘Vrnis only”. But here 
the error on E/A is rather large which was already fore¬ 
seeable from the also large error on energies of finite nuclei 
seen in[^ Besides E/A, the parametrization ‘Vrms only” 
yields also a very large error on the symmetry energy J. 
This happens because the selection of nuclei with radius 
data is smaller along isotopic chains thus carrying even 
less isovector information. 


Golumn 3 of figure shows the effect of variation of 
fit data on response observables in ^°®Pb. The three gi¬ 
ant resonances and the polarizability ao are known to 
have a one-to-one correspondence with each one NMP 
(see sections 6.3 6.5 and ref. [31]): the GMR with K, 
the GQR with m*/m, the GDR with ktrk, and ao with 
Osym- These pairs of highly correlated observables show 
exactly the same trends in the hgure which confirms the 
correspondence also from this point of view. Unlike NMP, 
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Fig. 8: Results on NMP (left two columns) and specific observables in finite nuclei (right two columns) for the three fits to 
different data sets: “SV-min” = fit to the full standard data pool of [^, “i? + rrms only” = energies and r.m.s. radii in the 
pool plus gaps and l*s splittings, “i? only” = energies in the pool plus gaps and l*s splittings. Left column: isoscalar NMP. 
Second from left colnmn: isovector NMP. Third column: Giant resonance energies and polarizability in ^°®Pb. Right column: 
binding energy and fission barrier in ^®®Hs and binding energy and a-decay energy in the the hypothetical super-heavy nucleus 
Z=120/N=182. 
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observables in finite nuclei allow a comparison with exper¬ 
imental data. GMR and ao are related to the static re¬ 
sponse properties K and J. Comparison with data shows 
clearly that inclusion of formfactor information drives into 
the right direction while a fit without i?diffr,c and ac yields 
less favorable values. This indicates that information from 
the formfactor carries correct physics about the nuclear 
response although these data were found in the previ¬ 
ous figure to stay slightly in conflict with energy informa¬ 
tion. The situation is mixed for the dynamical response 
properties. The GDR improves with omitting i?diffr,c and 
ac while GQR deteriorates this way. This is very likely 
pure chance because dynamic response is not embodied in 
ground state data. An explicit ht of dynamic response, as 
done e.g. in SV-bas [52] j is anyway recommended. 

The situation is mixed for the super-heavy elements 
(SHE) shown in column 4 of figure As several times 
before, all predictions agree with each other within their 
error bars. There remains to comment a few trends of 
average values. Skipping i?diffr,c and ac (parametrization 
“A-l-rrms only”) improves EbC^^^Hs) but spoils R/(^®®Hs). 
Skipping, further more rrms.c to parametrization “E only” 
does not make much of further changes. However, the 
parametrization only” spoils both i?B(^®®Hs) and 

i3/(^®®Hs). There are even less trends in the averages for 
the SHE Z=120/N=182, except for ‘Vi-mg only” which is 
obviously totally inappropriate for SHE. This is also sig¬ 
naled by huge error bars which “ri-ms only” produces in 
all four observables of SHE. The data basis of r.m.s. radii 
is obviously too weak to allow for far extrapolations. The 
examples demonstrate that statistical analysis combined 
with variation of model or conditions delivers very useful 
insight into the model. 


6.3. Systematic variation of nuclear matter prop¬ 
erties (NMP) 


In this section, we exemplify the strategy 3.2 (see section 
[^, namely systematic variation of a model parameter or 
property. We have argued in section 3.1 that the NMP as 
defined in table are fully equivalent to the parameters 
of the SHF functional and have the additional advantage 
that they carry an intuitive physical meaning. We follow 
here the strategy of [52] and define a base point with the 
following four NMP fixed: K = 234 MeV, m*jm = 0.9, 
J = 30 MeV, and ktrk = 0.4. These values were cho¬ 
sen because they lead to a good reproduction of the four 
basic response properties in ^°®Pb: GMR, GDR, GQR, 
and ao- A fit with these four NMP fixed delivers the 
parametrization SV-bas. Starting from this base point, 
we have produced four chains of parametrizations where 
one of the NMP was varied systematically while the other 
three were kept fixed. A detailed description of this set of 
parametrizations is given in [32] . In the following we use 
this set for the examples. 

Figure shows the trends of the four crucial response 
properties in ^°®Pb with the four NMP. The lower left 
panel shows the variation with K. Only the GMR moves 
with K while the other three response properties react 
only weakly. This makes it obvious that the GMR is 
uniquely related to K. The upper left panel demonstrates 
in the same way the close relation between GQR and 
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Fig. 9: Peak energies of giant resonances (GMR, GDR, GQR) 
and dipole polarizability q_d in for sets of parametriza¬ 

tions with systematically varied NMP from [32] as indicated 
on the a:-axes. 


m*/m. And similarly, we see the unique relation between 
GDR and ktrk in the upper right panel. The lower right 
panel shows the trends with J. Here it is the dipole po¬ 
larizability a£)(^°®Pb) which shows the strongest depen¬ 
dence. There is some trend for the GDR, but this remains 
smaller than the trend for ao versus J or for the GDR’s 
dependence on ktrk- Thus we conclude from these four 
variations that there is a one-to-one relation between the 
four response properties and the four varied NMP. These 
findings are conhrmed even more clearly by correlation 
analysis in section [675] and figure [T4| therein. 

One may wonder why a variation of L, the slope of sym¬ 
metry energy, is not considered. The point is that J and L 
are highly correlated with each other (covariance near 1) 
for the present SHF functional and pool of fit data. Thus 
a plot versus L shows the same trends as a plot versus 
J. This feature may change for extended SHF function¬ 
als which deliver more flexibility in the isovector part of 
density dependence. An example in that direction is the 
extension (35) which indeed reduces the correlation be¬ 
tween J and L as will be seen later in the middle panel of 
hgure 

Figure 10 shows the trends of fission barrier Bf and hs- 
sion lifetime Tf of ^®®Hs (for its definition see section 3.2) 
for varied NMP. The dominant influence comes here from 
m*/m. This is not surprising because hssion is driven by 
shell effects and these are highly sensitive to m*/m. The 
hgure is, however, to some extend misleading. It gives 
the impression that hssion is exclusively related to m*/m. 
This view is not conhrmed by correlation analysis mm 
where the covariance of Bf with m* jm is about 0.2 which 
is some correlation but not a dominant one. The reason 
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Fig. 10: Fission barrier -B/ in ^®®Hs (lower panels) and fission halflife r/ (upper panels) for sets of parametrizations with 
systematically varied NMP as indicated on the a;-axes. 


Z=50 isotopes N=82 isotones 



neutron number N proton number Z 

Fig. 11: Energy of the low-lying 2^ state in the isotopic chain 
of Sn (Z=50) and the isotonic chain with N=82. Shown are 
results for SkI3 | 1U4 | and a variation of pairing strength while 
keeping the SHF parameters of SkI3 fixed. Data taken from 

HD- 


been particularly fitted for that. Larger deviations occur 
near shell closures that is the region where collectivity is 
not so large such that the mapping of the collective path 
by a mere quadrupole constraint becomes insufficient HD- 
One should employ here variationally optimized paths as 
delivered by adiabatic TDHF, see e.g. [15]. The effect of 
varied pairing is dramatic showing that low-lying 2“'" states 
are a sensitive probe for the pairing functional. This is a 
feature which has not yet been fully exploited in calibra¬ 
tion strategies. The effect is particularly dramatic for the 
case of reduced pairing (purple line with filled triangles). 
The excitation energies shrinks occasionally to near zero. 
These are cases where pairing breaks down and it shows 
that the standard values of pairing strength are rather 
weak (i.e. close to the break-down regime). 


6.4. Impact of ab-initio data 


that there are several other agents which all have some 
impact on the fission path. The next larger influences 
come from the surface energy, spin-orbit term (also acting 
through shell effects), and pairing. The example demon¬ 
strates that systematic variation of force parameters is 
very instructive to reveal influences, but should be aug¬ 


mented by correlation analysis (section 6.5 and point 2.2 
in section^ to avoid premature conclusions. 

Figure ^1| shows the E 2 + energies of low-lying 
quadrupole states in the chain of even Sn isotopes and 
even N=82 isotones. We show here the effect of a variation 
of pairing strength on the results. Similar as for the fission 
path, the collective path for large amplitude quadrupole 
vibrations is most strongly influenced by pairing and grabs 
a piece of influence from many other ingredients. Thus we 
show here only the strongest effect and refer the reader 
to [41| for the many more subtle variations. First, we 
point out that SkI3 with standard pairing provides satis¬ 
fying results, particularly for Sn in the mid shell region. 
We see again that the SHF functional is capable to de¬ 
scribe at once also large amplitude motion without having 


The nuclear many-body problem is much more involved 
than the electronic one. The microscopic nucleon-nucleon 
interaction has a huge repulsive core at short distances and 
yet one cannot ignore long-range correlations from zero- 
sound modes. Moreover, it is not even correct to speak of a 
nucleon-nucleon interaction because nucleons are compos¬ 
ite particles from quarks and gluons, and these constitute 
a highly non-linear field theory. In spite of all these enor¬ 
mous complications, much progress has been made with 
nuclear ab-initio calculations during the last two decades 
miEolEI]. Time is coming close that we can refer to ab- 
initio calculations in bulk matter, perhaps also in finite 
nuclei m- Earlier attempts to use that for calibration 
of mean-field models were promising, but not yet fully sat¬ 
isfying, see e.g. the study in [39] on the basis of relativistic 
Brueckner-Hartree-Fock data of Ab-initio calcula¬ 

tions have still further improved since then. In particular, 
there is one point which seems to be well settled. This is 
neutron matter at very low densities. This stage is called 
a correlated Fermi gas (CFG). In the limit of low densities 
p, or Fermi momenta kp respectively, the CFG energy can 
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be parametrized in terms of the neutron kinetic energy as 


E 

A^neut 


(fc/) = C 


E 

N kin ,neut 


{kf) 


(33) 


where ^ is, in principle, function of kp. We are interested 
here in the limit kp —> 0 for which ^ becomes just one 
number. Many-body theory predicts a value ^(0) = 0.38 — 
-0.44 for the CFG [rOTlfT^ . 

It is straightforward to compute for the SHF func¬ 
tional. With the constituents of the total energy Q, we 
can deduce = (fkin + ■?Sk) /Po and thus 

e = 1 + . (34) 

^kin 

To estimate ^ in the low-density limit, we recall kp oc 
and discuss the limit in terms of p because the SHF 
functional ([^ is given that way. The kinetic density in 
homogeneous matter becomes fkin oc r oc p^/^ |109) . The 
leading term in the SHF functional is the two-body contact 
interaction cx p^. The density dependent term behaves 
as oc Successful parametrizations require a > 1. 

This term, therefore, vanishes faster than p^ in the limit 
p —>■ 0 and can be neglected. The kinetic interaction term 
behaves as pr oc p®/^ which, again, can be neglected. It 
remains ^ oc I-|-p^/p®/® —>■ 1 in the limit p —>■ 0. Thus 
we see in purely analytical manner that the standard SHF- 
functional cannot meet at all the requirement ^(0) « 
0.42. 

In order to allow for a description of ^ < I, we extend 
the SHF functional by adding to 5sk another density de¬ 
pendent term 


^Sk,dens3 


■^P-1/3 
*-^0 


Po + C 


p-1/3 
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Pi 


Po'/^ (35) 


E 

E 


U) 

E 


o 

0) 


15 

10 

5 

0 

-5 

-10 


> 

0) 


m 

lU 

o 

o 


0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

- 0.4 

- 0.6 

- 0.8 



D 

o 

0) 


40 

30 

20 

10 

0 

-10 

-20 


■D 

DC 

o 

0) 


0 

-10 

-20 

-30 

-40 

-50 

-60 



E 

I 

> 

(/) 



O 

u. 

O 

+ 


Fiff. 12: Average and r.m.s. errors on basic bulk properties, 
binding energy Eb, charge r.m.s. radius r^ms, charge diffraction 
radius and charge surface thickness (Tgurf for the series 

of fits with the additional density dependent term (351 and 
four NMP fixed. Averages are taken over the pool of fit nuclei. 
The average deviation from data is shown by filled boxes and 
the r.m.s. deviation by the error bars around the boxes. The 
series is: “SV-min” = fit of standard functional to full data 
pool from 


-fCy = fit as SV-min but including the 


new term (351 in the functional; “-I-CFG” = fit as “-(-Cy 
but including the CFG limit as data point. 


It has the form of the density dependent term in the func¬ 
tional § with a = —1/3. Handled as additional term, 
it allows to maintain the quality of the standard ansatz 
while opening the chance to accommodate the CFG limit. 

We have produced two new parametrizations with the 
extended functional, now including also the term (351. 
Starting point is here, again, SV-min, a fit to the standard 
pool of fit data. The first new set is “-(-(7^ and uses 
the same data as SV-min. The second new set “-I-CFG” 
adds to these data the CFG condition ^(0) = 0.42 as fur¬ 
ther data point. 

Figure shows the resulting average and r.m.s. er¬ 
rors for the basic blocks of observables. The averages stay 
at zero level throughout for Ep and vary somewhat for 


radii and surface thickness. The fit with more 

freedom in density dependence, yet without enforcing the 
CFG point, leads always to lower average and r.m.s. er¬ 
rors. This confirms the impression gained in previous sec¬ 
tion that the slight conflict between energy and formfactor 
calls for improved density dependence. However, the step 
to “-I-GFG” deteriorates again radii and surface thickness 
in the average approximately back to the status of SV- 
min. The freedom coming along with the new term has 
been fully exploited by the GFG point. 

Figure shows NMP and response properties in 
208pb fQj. parametrizations introduced in this section. 
There are marked changes in the results from the three 


parametrizations. But they are all compatible within the 
error bars. Nonetheless, it may be instructive to discuss 
the changes in the center values. Similar to the previous 
figure, we find that the results of SV-min and “-I-GFG” 
are closer together while makes generally the 

largest changes. Exceptions are the pair kqdr with GDR 
and K with GMR where “-I-CFG” produces too low val¬ 
ues. The strong impact of the GFG point is plausible for 
K and GMR because these observables have direct rela¬ 
tion to density dependence. The fact that the GFG point 
pulls the GMR away from the data indicates that the ex¬ 
tended model is not yet ideal, although the data may be 
corrected by fitting with constraint on K in the manner 
as was done for SV-bas [32]. It is also noteworthy that 
the error bars for the free fit of the extended density de¬ 
pendence in “+Cp produces particularly large un¬ 

certainties for K and GMR, the two observables sensitive 
to density dependence. This large uncertainty is much re¬ 
duced by the step to “-I-GFG” which indicates once more 
that the GFG point has a strong impact on the density 
dependence of the model. 

The example is also quite revealing as it shows that 
some long-range correlations (at least those at extremely 
low density) are not incorporated into the SHF functional. 
This is annoying because these corralations are smooth 
with mass number and thus should be contained in a 
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Fig. 13: Results on NMP (left two columns) and response prop¬ 
erties of (right column) for the series of three forces as 

in figure 


smooth functional. Again, we see that we need to improve 
on the density dependence of the SHF functional. 


6.5. Correlation analysis 

In sectioncorrelations, or covariances respectively, be¬ 


tween observables had been defined with eq. (32d). In this 


section we will give a few examples for a small, but rel¬ 
evant, selection of observables. The selection covers: the 
four NMP K, m*Im, J, ktrk together with the four re¬ 
lated response properties in ^°®Pb, namely GMR, GDR, 
GQR, and an] for static isovector response we look also 
at L, neutron skin rskin(^°®Pb), and slope of the neu¬ 


tron equation of state dpE/Aneut 


p=0.1 


; finally we check 


extrapolations to exotic nuclei, Eb for the very neutron 
rich ^'^°Sn and Eb as well as Qa for the super-heavy ele¬ 
ment Z=120/N=182. Three parametrizations will be con¬ 
sidered: SV-min as the standard, “E only” from section 
6.2 to exemplify the imp act o f the choice of fit data, and 
“-t-Cp from section 6.4 to explore the impact of an 

extended functional. 

Figure shows the correlations between the selection 
of observables as indicated on the axes. As every observ¬ 
able can be related to any other, we show it as correlation 
matrix. Let us start with the parametrization SV-min, the 


SV-min 


Q„(120/182) 

Eb(120/182) 

Egf'^Sn) 

n.skin(^“Pb) 

aD(^“Pb) 
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Fig. 14: Correlation matrix for a selection of observables as in¬ 
dicated. Three cases are shown: SV-min = fit of standard SHF 


functional to the standard pool of data from |32| . “-|-Gy 


-l/3„ ^ 

unconstrained fit as SV-min but with the functional extended 
by the density dependent term (351, “E only = fit of stan¬ 


dard functional to standard data pool without radii and surface 
thickness (see section [G^. 


21 




























































































straightforward fit to the standard pool of data [32]. We 
see clearly the strong correlation between K and GMR 
(static isoscalar response), m*/m and GQR (dynamic 
isoscalar response), ktrk and GDR (dynamic isovector 
response), and for the whole group of static isovector re¬ 


sponse covering here the NMP J, L, d/dpE/A 


neut 


(see 


section [3^ together with the finite nuclei observables ajj 
and Tskin (labeled “n.skin” in the figure). For the exotic 
nuclei, there is some correlation of ^^^(^^^Sn) to isovec¬ 
tor static response which is plausible for a extremely neu¬ 
tron rich nucleus with very large isospin. The two data 
in the super-heavy element Z=120/N=182 are nearly de¬ 
coupled from the other groups of data. One may admit 
some correlation with iF&GMR. This super-heavy region 
is dominated by shell effects and not so much by the bulk 
properties included in this correlation matrix. It is also 
interesting to remark that Eg and Qa are almost uncor¬ 
related. Differences of energies often filter quite different 
aspects of a model than the energy as such. 

We now look at the correlations produced by the next 
two parametrizations, with less data in “E only” and with 
more terms in the model in . In both cases, 

we see generally reduced correlations what we could have 
expected. It is, however, remarkable how small the re¬ 
duction is leaving the gross structure of the correlations 
fully alive. There remain clearly the blocks for the four 
different responses. Larger changes are seen for isovector 
static block. And here, the changes are very different for 
the two different variations of the fit. For example, while 
“i? only” degrades the correlations for J thus making L 
the leading bulk parameter for isovector static response, 
the other parametrization degrades just L and 

leaves J more intact. In such details, we see indeed a de¬ 
pendence on model and data pool. But recall that the 
general structure of correlations was robust. 

There is another interesting detail concerning the super¬ 
heavy element Z=120/N=I82. The independence of Eb 
and Qa persists. But there is now more correlation be¬ 
tween Qa and AT&GMR than was seen for SV-min. This 
shows that correlations are not necessarily reduced by less 
data or a large model. It can also happen that more free¬ 
dom allows to explore connections which were hindered 
before by a too rigid model. 


6.6. Spin sensitive applications: spin modes, odd 
nnclei, TDHF 

This section is devoted to the spin properties of SHF. 
Spin properties are crucial for a broad range of applica¬ 
tions. They play a role already for ground state properties, 
namely in case of odd nuclei, see e.g. [TT0l[Tm[TT2l[TT3] . 
They are key players in spin excitation modes, see e.g. 
fTTTlfnmiTTCI . they play a role in rotating nuclei, see e.g. 
[miinsi, and they can contribute to energetic heavy- 
ion collision described by time-dependent Hartree-Fock 
(TDHF) with the SHF functional, for an example where 
the impact of a spin excitation is discussed explicitly see 
m- Unfortunately, the spin terms are presently the least 
well controlled part of the SHF functional. Their proper 
calibration is yet an open problem and they are most prone 
to introduce instabilities into the SHF equations, for a crit¬ 
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Table 3: Restrictions on the parameters of the general 
SHF functional ([^ in case that the functional is derived 
from the Skyrme force ansatz (36). 


ical analysis see [251 (23 [2H1 US] ■ The development of the 
spin branch in SHF is still in its first stages, this means 
in the language of section [5] at the preparatory level of 
point although one can find also some steps of the error 
analysis according to point and sub-points |3.I| together 
with |3.4| therein. This latter aspect will shine up briefly 
in the subsequent example of odd nuclei and of TDHF. 
But generally, the aim of the present section is mainly to 
demonstrate the capabilities of SHF to deal also with spin 
properties. 

Before going on, we have to look at the spin aspects of 
the basic SHF functional ([^, so far not discussed. Some 
spin terms, namely the time odd spin-orbit terms oc 
in the time-odd part of the functional (6c I, are fixed by the 
requirement of Galilean invariance. Yet open remain the 
parameters of the genuine spin terms and the tensor 
spin-orbit terms oc C^. From a strictly density functional 
perspective, one can consider these all as free parameters 
of the model. This leaves still a large number of 7 more free 
parameters whose independent calibration would require 
access to a great amount of reliable, spin-sensitive data. 
Alternatively, one can approach the problem from a formal 
side. The density-matrix expansion sketched in section [23] 
yields not directly an energy functional, but rather a zero- 
range effective interaction, the much celebrated Skyrme 
force. For its detailed form and a critical discussion of the 
notion of an effective interaction see appendix]^ Deriving 
the SHF functional from the Skyrme force yields a well de¬ 
fined connection between the spin terms and the leading 
terms, see table]^ These restrictions on the spin-orbit pa¬ 
rameters from the force definition of the SHF functional 
are often discarded. Many parametrizations simply ignore 
the tensor spin-orbit terms, thus supposing Cq = 0 and 
Cf = 0. There are also reasons to override the restric¬ 
tion on the spin-orbit parameter CY'^ = |I04j . All 

functionals used in this paper ignore the tensor spin-orbit 
terms and treat CY"^ as independent parameter. The spin 
gradient terms oc when taken as given in table lead 
for most functionals to dramatic instabilities. They must 
be discarded. What remains then from the force concept 
is the parameter connection from Galilean invariance and 
the determination of the pure spin terms cx C^- Galilean 
invariance is the unavoidable minimum requirement and 
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Fig. 15: Difference of theoretical and experimental binding 
energies (lower panel) and neutron odd-even staggering A® 
(upper panel) along the isotopic chain of Sn nuclei computed 
with the parametrization SV-bas m- The -Emin stands for the 
minimal Galileo invariant functional and Ej^j^ ^2 for the the 
functional including additionally the terms according to the 
force definition, see table Data extracted from m- 


will always be used. Above that we have the freedom to 
employ the spin and/or tensor spin-orbit terms. Although 
one cannot take the full “force” terms consistently into 
account, it is convenient to exploit it for the spin terms 
and the tensor spin-orbit part yielding the connections as 
given in table[^ As this is an ad-hoc decision, justification 
can come only from experience. 

As a first example, we address odd nuclei picking a key 
result from |113j . Figure [I^ shows the deviation of binding 
energies from experimental values and the odd-even stag¬ 
gering (191 along the chain of Sn isotopes. Two stages of 
the functional are considered, the minimal Galilean invari¬ 
ant form Emin corresponding to = 0 and the inclusion 
is the cr^ terms with parameters as given in table The 
trend of the binding (lower panel) is generally very similar 
for odd and even isotopes. There is a good reproduction 
in the region N = 70 — 82. The larger deviations mid 
shell indicate that these nuclei require collective correla¬ 
tions (not included here) for a proper description. Odd 
nuclei make no difference in that respect. Of course, both 
variants yield the same values for even isotopes because 
the time-odd terms are inactive there. It remains to look 
at differences for odd isotopes. The effect of the spin terms 
is small, of the order of 0.2 MeV and in the majority of 
cases enhancing slightly the deviations. Thus the results 
for binding energies give a slight preference for the Emin 



CO [MeV] 


Fig. 16: Isoscalar Ml strength distribution for the spin-flip 
mode in deformed ^®®Gd and the spherical ^°*Pb computed 
with separable RPA using the parametrization SV-bas. The 
RPA results are shown by heavy line and the strength dis¬ 
tribution from pure two-quasiparticle excitations (from static 
mean field only) by faint lines. The experimental distribution is 
shown for ^®®Gd. The remarked experimental resonance peak 
is indicated by an arrow for ^°®Pb. Adapted from [HQ]. 


model. 

The upper panel of figure shows the odd-even stag¬ 
gering. It is very satisfying that the parametrization SV- 
bas provides at once a satisfying reproduction of the data. 
At a more detailed level, we see differences between the 
two choices for the spin terms. There are cases where 
£'min,o -2 performs better (particularly at the side of small 
N). But, again, we find that the majority of isotopes 
looks better with the minimal Galilean invariant choice 
This slight preference for Emin was also found for 
other examples and other parametrizations |113j . But the 
differences are probably too small to be decisive and more 
studies are still necessary to settle the case. Recall that 
the most important result of this brief demonstration is 
that the SHF functional yields immediately an acceptable 
description also of odd nuclei. 

Figure [TH] shows as next example a result for spin exci¬ 
tations, here the Ml strength distribution of the spin-flip 
mode, computed including the spin terms of te SHF func¬ 
tionals. Two rather different nuclei are considered, the 
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well deformed mid-shell nucleus ^^®Gd and the spherical, 
doubly magic The RPA result is shown together 

with the mere two-quasiparticle excitations (i.e. without 
the dynamical effects from the RPA residual interaction 
Bnn). One sees a marked upshift of strength by RPA (also 
called collective shift) due to the repulsive residual inter¬ 
action in the isoscalar spin channel. This upshift brings 
indeed the strengths into the right energy range. Thus 
far we see once more that a proper SHF functional is also 
appropriate for spin modes, at least what gross features is 
concerned. Looking closer at the case, we see a problem. 
The experimental distribution differs between deformed 
and spherical nuclei. It has two pronounced peaks for 
the deformed ^®®Gd but only one for the spherical ^°®Pb. 
The theoretical result shows a two-peak structure for both 
cases. This has been studied extensively in |116LI120] for a 
broad selection of SHF parametrizations. There are some 
parametrizations which produce one single peak for ^*^®Pb, 
but then one peak only also for ^®®Gd. With presently 
available functionals it seems impossible to describe both 
nuclei equally well. The source of the problem lies most 
probably in an incomplete inclusion of spin-orbit terms. 
A satisfying solution has yet to be worked out. 

Finally we demonstrate the enormous capabilities of 
self-consistent mean field theories by examples from time- 
dependent Hartree-Fock (TDHF). TDHF allows to sim¬ 
ulate a wide range of nuclear processes covering reso¬ 
nance dynamics, non-linear excitations, and particularly 
large amplitude motion with substantial rearrangement of 
structure in the course of time as it occurs, e.g., in heavy- 
ion reactions. First nuclear TDHF calculations came up 
already shortly after the appearance of static SHF cal¬ 
culations m- Two decades ago, TDHF has been re¬ 
vived because computing capabilities had evolved suffi¬ 
ciently far to allow realistic large-scale calculations, for 
examples of a couple of very different applications see 
[1221112311124111^5] . A combination of TDHF with density- 
constrained Hartree-Fock even allows to compute sub¬ 
barrier fusion cross sections [TMllT?f| - Stepping forward 
to still larger scales, TDHF is now being applied to sim¬ 
ulate matter under astro-physical conditions in neutron 
stars and super-nova explosions For recent re¬ 

views on TDHF with SHF functionals see [13011131] . In all 
these applications, one found that SHF functionals provide 
at once a fair description of the processes studied. There 
are open ends though because spin terms can play a role 
in energetic collisions and we have seen above that just 
these terms are not yet well settled. We will exemplify in 
the following both, the capabilities and the uncertainties 
of TDHF simulations. 

Figure [TT] shows snapshots of the Wigner function from 
a TDHF simulation of a central collision ^®0-|-^®0 [132j . 
The simulation was performed on a Cartesian 3D mesh, 
for details of the physics and numerics of TDHF and 
the 3D code see m- The Wigner function is con¬ 
structed from the one-body density matrix g(r, r') by 
Fourier transforming with respect to y = r — r', i.e. 
f{r,k) = J -1- yl2,r — y/2). Originally in¬ 

troduced in [133| . it is the attempt to provide a phase- 





t = 67.8fm/c 



z[fml 

t = 155.8fm/c 



zffml 


space picture of a quantum state, for details see [134| . The Pw- If- Snapshots of the Wigner distribution in the z-fcj-plane 
Wigner function has the weakness that it is not positive hom a TDHF simulation of a central ^0+ ®0 collision using 
semi-definite, thus preventing a strict probabilistic inter- parametrization SkI3 Adapted from 
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Fig. 18: Time evolution of the quadrupole momentum Q 20 for 
a central collision '*‘^Ca+'^°Ca computed for the parametriza- 
tion SLy4 with three different options for the spin terms 
as indicated. Adapted from | 138 |. 

pretation. Nonetheless, it can serve as a useful illustration 
of a dynamical process and the availability of all detailed 
wavefunctions in a TDHF calculation allows us to produce 
this insight. The Wigner function in a 3D calculation is 
six-dimensional and thus rather difficult to handle. We 
thus look at a 2D cut through the 6D function, namely 
along the reaction axis z and associated wavenumber k^. 
The snapshots are sorted from above to below in order 
of increasing time. The initial stage shows nicely the two 
Fermi spheres of the nuclei, unperturbed except for 
the c.m. boost to positive kz for the left nucleus and to 
negative kz for the right one. The second panel shows the 
situation at the point of closest contact. Although the cor¬ 
responding local density represents a compact compound 
system, the phase space picture shows that the compound 
stage has not lost its memory at the initial state. There 
are two distinct fragments in phase space. This feature 
persists to the next plot the fragments separating. They 
are strongly perturbed, but kept their identity from the 
initial nuclei. And in the final stage (lowest panel), two 
fragments are flying apart, internally excited and with re¬ 
duced c.m. motion, but still keeping the two-fragment 
structure of the initial state. This indicates that pure 
mean-field motion does not allow full melting into one 
thermal compound state. This becomes possible only if 
two-body collisions are added to the description. This is 
achieved in the semi-classical framework by the Vlasov- 
Uehling-Uhlenbeck scheme |135) . Adding a collision term 
to quantum mechanical TDHF turns out to be much more 
involved, although just recently promising attempts have 
come up m- This example indicates that TDHF dy¬ 
namics is likely to underestimate dissipation processes. It 
still provides relevant guidelines for the gross structure of 
a process and can be even quantitatively correct in the en¬ 
trance state of a heavy-ion reaction where nucleon-nucleon 
collisions did not yet have a chance to play a role. For 
example, one can describe very well fusion cross sections 
with TDHF, see e.g. m- 

Figure [TS] shows the result from a TDHF calculation for 
a central collision of ^°Ca-|-^°Ca at relative center of mass 


energy 240 MeV |138j . Shown is the quadrupole momen¬ 
tum along the reaction axis as function of time. This quan¬ 
tifies the relative distance of the fragments and illustrates 
the type of reaction as well as damping of c.m. motion. In 
this case, three different choices for the SHF functional are 
compared: The parametrization SLy4 just with the min¬ 
imal Galilean invariant terms, SLy4 with the cr^ terms 
included, and SLy4 with cr^ plus tensor spin-orbit terms 
included. Unlike the case of ground states of odd nuclei, 
we see in this highly dynamical case a great sensitivity to 
the spin and spin-orbit terms. Their inclusion turns a fu¬ 
sion dynamics (heavy line for SLy4 without spin terms) to 
inelastic scattering which means that they reduce dissipa¬ 
tion. It ought to be mentioned, however, that the leading 
spin orbit term cx acts in different direction, signif¬ 
icantly enhancing dissipation |1391 II40] . Thus there is a 
very subtle interplay of the various terms in the SHF func¬ 
tional. With the new generation of computers and codes, 
one will now be in a position to systematically investigate 
the intrinsic structure of heavy-ion reactions. And, as in 
all previous example of this section, we emphasize that 
the SHF functional delivers in general a relevant picture 
of structure and dynamics also for cases where spin plays 
a role. This is a safe basis for then going into details. 

7. Concluding remarks 

The Skyrme-Hartree-Fock (SHF) model for a self- 
consistent mean-field description of nuclear structure and 
dynamics has been inspected from different perspectives: 
formal background, range of applicability, performance, 
estimates of extrapolation uncertainties, and correlation 
(covariance) between different observables. SHF is based 
on an effective energy-density functional which is formally 
motivated by a low-momentum expansion of a most gen¬ 
eral effective interaction and its free parameters are cal¬ 
ibrated by a fit to empirical data, mostly from nuclear 
ground states. This mixed origin calls for particularly 
careful exploration of the capability, reliability, and limits 
of the model. One key task is to estimate the uncertainties 
in extrapolations. Related to that is the other task to find 
out the dependencies between observables (and model pa¬ 
rameters) which helps to identify, e.g., bottlenecks in the 
model and to cure them. A great manifold of strategies to 
attack these task was outlined in section A large part 
of them was exemplified in the subsequent result section. 
The same pool of fit data was used in all examples to ren¬ 
der the tests consistent with each other and so maintain 
comparability. These data were bulk properties of nuclear 
ground states in chains of semi-magic nuclei selected to be 
good mean-field nuclei (only very small effects from collec¬ 
tive ground-state correlations). We summarize the basic 
findings of that survey in key sentences: 

1. There is no single ideal strategy to estimate the un¬ 
certainties in predictions from a model. A bundle 
of different methods has to be invoked to develop a 
sound intuition for the limits and capabilities of a 
model. Statistical analysis in connection with fits 
is one important part. Harder is an estimate of the 
systematic error. It has to be approached from differ¬ 
ent sides, formal analysis as well as variations of the 
model, and yet will always leave open ends. 
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2. Concerning the basic performance, modern SHF func¬ 

tionals have reached a high quality of description of 
ground state properties in the whole known nuclear 
landscape, from light nuclei, as up to super¬ 

heavy elements. For example, energies are reproduced 
with a r.m.s error of about 0.6 MeV. 

3. Although not explicitly fitted for that, modern SHF 
functionals can be used successfully for many other 
observables and processes as resonance excitations, 
astro-physical nuclear matter, nuclear rotation at 
high spin, and large-amplitude collective motion (low- 
energy quadrupole vibrations, Hssion, fusion, heavy- 
ion scattering). The latter example goes, in fact, often 
beyond a mere mean-field treatment. 

4. Although small, the residual errors are not statisti¬ 
cally distributed, but show still unresolved trends as 
typical fluctuations between doubly-magic and mid¬ 
shell nuclei and, more puzzling, a small but sys¬ 
tematic increase of under-binding towards heavy and 
super-heavy nuclei. The latter trend is probably due 
to yet insufficient modeling of density dependence of 
the SHF functional. 

5. Extrapolation errors from statistical analysis provide 
a lower estimate (yet without contributions from sys¬ 
tematic errors) which is already extremely useful to 
distinguish safe and risky extrapolations. Not sur¬ 
prisingly, errors increase with increasing distance to 
the 6t data. For example, extrapolations to super¬ 
heavy elements are still found to be reliable while 
error bars grow huge when predicting properties of 
neutron stars. On the other hand, interpolations in 
the vicinity of the fit data are safe and yield devia¬ 
tions as small as the average errors in the ht pool. 

6. There is a hierarchy of importance within the terms of 
the SHF functional. From a formal side, importance 
shrinks with increasing order of derivatives involved. 
From the empirical side, isoscalar properties are bet¬ 
ter determined than isovector ones due to the lim¬ 
ited extension of known nuclei along isotopic chains. 
Taking the most important terms delivers a minimal 
model with only 6 free parameters which provides al¬ 
ready the high quality of ground state properties as 
the final model employing all 11 free parameters of the 
SHF functional. The further terms above the mini¬ 
mal model serve to optimize further features as, e.g., 
giant resonance excitations. 

7. The impact of the fit data was checked by perform¬ 
ing fits with successively omitting a group of data. 
From the four groups of fit data (binding energy Eb, 
charge r.m.s. radius r^ns) charge diffraction radius 

charge surface thickness ac), Eb is found to 
be the most decisive one. Fits with Eb alone yield 
already an acceptable model. Adding Tmis to the pool 
of Eb comes another step closer to the full ht while 
considering Crms alone is insufficient. The small drifts 
of average errors and of predicted observables when 
adding i?diffr and ac shows that each group pulls the 
ht into a slightly diherent direction. This conhict be¬ 
tween observables points to a rigidity in modeling the 


nuclear surface prohle which, again, is probably re¬ 
lated to the modeling of density dependence. 

8. A particularly instructive analyzing tool is the sys¬ 
tematic variation of model features, here exemplihed 
by the variation of nuclear matter parameters (NMP). 
It illustrates the impact of such a parameter on an 
observable of interest. This allowed to demonstrate 
the near one-to-one correspondence of NMP and giant 
resonances in ^°®Pb, namely incompressibility K with 
the giant monopole resonance, isoscalar ehective mass 
m*/m with the giant quadrupole resonance, symme¬ 
try energy J with the dipole polarizability ao, and 
the TRK sum-rule enhancement ktrk (alias isovec¬ 
tor ehective mass) with the giant dipole resonance. 
A much diherent behaviour is seen for hssion barrier 
and lifetime in the super-heavy nucleus ^®®Hs. These 
observables gather inhuences from many diherent pa¬ 
rameters and have no prevailing correlation. 

9. Ab-initio data are not yet regularly included in the 
calibration of the SHF functional. As an example for 
ab-initio input, the neutron gas at very low density 
was considered. This data point can meanwhile be 
computed very well with the theory of the correlated 
Fermi gas (CFG) and it is exclusively determined by 
long range correlations from zero-sound modes of the 
gas. It turns out that the standard SHF functional is 
not able to accommodate this data point from CFG. 
A new density dependent term must be added to allow 
an adjustment which then works acceptably well. The 
extended model allows also to test the ehect of an ad¬ 
ditional density dependence when htting to the stan¬ 
dard pool of data (excluding the CFG point). This 
ht resolves the conhict between observables discussed 
under point thus conhrming that modeling of den¬ 
sity dependence is still an issue. 

10. Correlation (or covariance) analysis was exemplihed 
for a small selection of relevant observables cover¬ 
ing the four essential groups of response properties, 
isoscalar static, isoscalar dynamic, isovector static, 
and isovector dynamic. This conhrmed very clearly 
the one-to-one correspondence between NMP and gi¬ 
ant resonance, already worked out under point 
Noteworthy is a large block of isovector static observ¬ 
ables covering, beside J and au, the neutron skin, 
pure neutron matter, and the slope of symmetry en¬ 
ergy L. All these relations are rather robust under 
changing conditions of modeling and selection of ht 
data. 

11. Spin terms in the SHF functional leave still open ques¬ 
tions. A part of them is determined by the require¬ 
ment of Galilean invariance. Other terms can be hxed 
by the concept of a Skyrme “force”. But this set¬ 
ting often leads to principle instabilities in the spin 
channel. Particularly dangerous are the spin-gradient 
terms which thus are discarded in all applications. 
There remains, again, a decision based on empirical 
data. The case is not yet settled and requires further 
research. 
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12. The pairing part of the functional has not been dis¬ 
cussed in this paper. Similar as the spin terms, it 
is also not yet fully explored. The present form is 
optimized to reproduce the odd-even staggering in 
isotopic and isotonic chains of semi-magic nuclei and 
serves to produce well defined ground states for open 
shell nuclei. Pairing is found to have a large impact 
on low-lying quadrupole excitations and fission. This 
may require still a further fine tuning for such pro¬ 
cesses. 


13. As mentioned in point the SHF functional is of¬ 
ten successful even in applications beyond pure mean 
field. Such methods require the energy overlap be¬ 
tween different mean field states. These can be eval¬ 
uated with the SHF functional as long as the differ¬ 
ent states remain close (Gaussian overlap approxima¬ 
tion). However, projection onto good particle number 
or angular momentum require also overlaps between 
remote states and the SHF functional in its standard 
form become inapplicable for this task. New model¬ 
ing based strictly on a force concept is needed, and 
underway to resolve this problem. 


Altogether, the SHF functional has proven to deliver an 
extremely useful and reliable reproduction of static and 
dynamic features of nuclei as far as they are accessible 
to a mean-field description. One could be highly satis¬ 
fied with what has been achieved. However, the success 
shows the capability of SHF and motivates to ask for more. 
Such more is also highly desirable in view of the the re¬ 
markable progress in producing exotic nuclei and of the 
far extrapolations implied in astro-physical applications 
(extremely neutron rich isotopes in the r-process, neutron 
stars, nuclear matter under supernova conditions). Work 
on further scrutinizing and improving the SHF functional 
along the lines summarized above is in progress in many 
research groups around the world. 
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A. The Skyrme force 


The SHF functional was motivated in section [2)3| through 
the density-matrix expansion gniisniisi]- This expansion 
acts on a supposed microscopic interaction, the T ma¬ 
trix from some ab-initio model, and delivers the “Skyrme 


force” in the form (tensor terms are ignored here) 

Vsk = to(l+2;o-P<7)<5(ri2) 

+ %{^+X3p„)p°‘ (ri)(5(ri2) 

-f^(l-fcciP^) (^6{ri2)k'^ + fc'^(5(ri2)) 

-|-t2 ( 1 +2^2 At ) (r 12) fe 

-|-it4(d-i -1- 0 - 2 ) • fc' X i5(ri2)k , (36) 

ri2 = ri - r2 , 

P<T = ^{^ + ^ 1 ^ 2 ) , 

fe = -^ (Vi - V 2 ) , fe' = ^ (Vi - V 2 ) . 


where k acts to the right and k' to the left. We put the 
notion “force” in quotation mark because this object de¬ 
pends on the density which is produced by the wave func¬ 
tion on which this force acts. This is not a standard two- 
body interaction, but an effective force designed for the 
only purpose to derive an energy functional (for a detailed 
discussion of the principle problems see [69] 1. The eval¬ 
uation of the expectation value of Vsk with a BCS state 
yields indeed a functional of the form (§. In fact, it yields 
more than that. The 10 parameters of the “force” (36) are 
distributed over the 23 parameters of the SHF functional 
(when counting each term with a separate free parameter). 
This establishes a couple of relations between the param¬ 
eters of the SHF functional. They embrace, of course, the 
restriction from Galilean invariance already embodied in 
the time-odd part (6c) and add seven more relations which 
fix all spin terms as shown in table 

It is interesting to note that the Skyrme force is par¬ 
ticularly restrictive with respect to the spin-orbit terms. 
The tensor spin-orbit terms oc are fully determined by 
the kinetic terms as is given by Cf and and the 
CY'^ = CY^ ties the isovector to the isoscalar spin-orbit 
term. The first restriction is neither harmful nor beneficial 
as this tensor spin-orbit term can easily be compensated 
by the normal spin-orbit term. The second restriction 
is more serious and to some extend questioned by phe¬ 
nomenology nng. Most traditional parametrizations em¬ 
ploy this restriction. Recent fits, particularly those follow¬ 
ing [32] discard the restriction to maintain independence 
of the isovector spin-orbit parameter. 


B. Pool of fit data 


Table summarizes the pool of data as defined in [32] ■ 
These data were used to determine SV-min, SV-bas and 
the series of forces with systematically varied NMP (used 
here, e.g., in section 6.3). They are used for the new fits 

The detailed 


in this paper, see e.g. sections |6.2| or |6.4| 
values for the listed observables are found in the original 
paper [53] . 
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